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COS  ATI  COOES  1 

FIELD 

GROUP 

SUB-GROUP  1 

Several  tasks  of  research  on  turbulent  flames  have  been  pursued.  The  focus  has 
been  on  low  to  moderate  Reynolds  number,  premixed  flames  which  are  dominated  by 
flamelet  structures,  i.  e.,  the  major  heat  releasing  chemical  reactions  occur  in  thin  sheets. 
Large  density  changes  occur  across  flamelets  and  therefore  it  is  expected  that  these 
structures  have  an  important  influence  on  turbulence  generation.  The  focus  of  the  research 
has  thus  been  on  flamelet  characteristics. 
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19.  Abstract  (continued) 


A  fractal  representation  of  surfaces  in  turbulent  flows  is  found  to  be  a  very  useful 
tool  for  modeling  purposes.  Experiments  support  the  hypothesis  that  flamelet  surfaces  can 
be  represented  by  fractal  surfaces,  but  more  experiments  are  need  before  the  evidence  can 
be  considered  conclusive.  Passage  time  pdfs  which  reflect  the  dynamic  character  of 
flamelets  are  also  useful  in  combustion  modeling.  Such  distributions  were  measured  and 
found  to  be  reasonably  well  represented  by  one  of  two  proposed  functional  forms  for  the 
pdfs.  If  a  very  exact  representation  is  required  neither  function  seems  ideal.  Flamelet 
position  and  conditioned  velocity  measurements  were  performed.  The  pdf  of  flamelet 
position  is  nearly  guassian  and  the  standard  deviation  of  the  distribution  is  a  measure  of 
flame  brush  thickness.  The  velocity  measurements  reveal  the  flamelet  to  be  a  source  of 
flow  acceleration.  But  experimental  problems  have  hindered  the  gathering  of  quantitative 
data  on  the  magnitude  of  this  acceleration. 
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I.  INTRODUCTION 

The  original  objective  of  this  work  was  to  employ  laser  velocimetry  (LV),  laser 
induced  Rayleigh  scattering  (RS)  and  pressure  probe  techniques  to  measure  velocity,  gas 
density  profiles  and  pressure  in  premixed  turbulent  flames  for  determining  the  influence  of 
combustion  and  density  fluctuations  on  turbulence.  The  purpose  was  to  study  the 
generation  of  turbulence  and  the  enhancement  of  turbulent  transport  in  the  reaction  zone  of 
turbulent  flames  in  order  to  improve  our  understanding  of  these  processes  and  to  improve 
turbulent  combustion  models.  A  primary  focus  of  the  experiments  was  to  measure  terms 
appearing  in  the  turbulence  energy  balance  equation.  Two  terms  were  of  particular  interest 
--  the  term  representing  work  by  the  turbulent  shear  stresses  against  the  mean  flow  and  the 
term  giving  turbulence  generation  by  the  mean  pressure  gradient 

The  experiments  were  to  be  carried  out  on  a  laboratory  v-flame  burner  which  had 
been  studied  previously  at  Cornell  [1,2].  The  burner  supports  flames  held  by  rods  in  grid 
turbulence.  With  this  burner,  low  turbulence  Reynolds  numbers  are  obtained  [3],  and 
reaction  is  confined  to  thin  sheets  for  the  conditions  studied  [2].  The  burner  is  described  in 
more  detail  below. 

As  the  research  progressed  there  were  changes  in  the  objectives  as  a  result  of  new 
findings  both  at  Cornell  and  elsewhere,  new  capabilities,  and  of  new  ideas. 

During  the  period  of  funding  the  following  research  was  accomplished,  funded 
solely  or  partially  by  ARO:  1 )  development  of  models  for  premixed  turbulent  flames  based 
on  a  fractal  description  of  reaction  sheets,  2)  extension  of  fractal  based  modeling  concepts 
to  jet  mixing  and  turbulent  jet  diffusion  flames,  3)  measurement  of  passage  times  in 
premixed  turbulent  flames  using  Rayleigh  scattering  with  results  compared  to  the  model 
proposed  by  Bray,  Libby  &  Moss  [4],  4)  conditional  LV  velocity  measurements  in  v- 
flames,  and  5)  preliminary  studies  of  the  fractal  character  of  v-flames.  Much  of  this  work 
has  been  reported  previously  in  journal  articles  [3,5]  or  Cornell  University  reports  [6,7]. 

In  this  report  our  research  findings  will  be  summarized  and  discussed.  Emphasis  will  be 
placed  on  results  not  previously  reported.  It  should  be  noted  that  work  on  conditional 
velocity  measurements  continues  and  the  student  working  on  the  project  expects  to 
complete  his  degree  work  by  August  1987.  Modeling  work  also  continues  and  a 
manuscript  for  journal  submittal  is  in  draft  stage. 


II.  MODEL  DEVELOPMENT 


Models  for  predicting  the  turbulent  burning  velocity  and  the  mean  fuel  consumption 
rate  per  unit  volume  of  premixed,  turbulent  flames  have  been  proposed.  Fractal  ideas  have 
been  applied  to  the  modeling  of  jet  mixing  and  of  jet  flame  length.  This  work  is 
summarized  in  this  section. 

IL1  Model  for  Turbulent  Burning  Velocity 

For  a  range  of  turbulence  conditions  including  conditions  of  practical  interest, 
reaction  in  premixed  turbulent  flames  occurs  in  thin  reaction  sheets  referred  to  as  flamelets. 
For  modest  levels  of  turbulence  the  structure  of  these  sheets  approaches  that  of  steady 
laminar  flames.  For  this  condition  Damkohler  has  suggested  that 


(1) 


where  ut  and  u#  are  the  turbulent  and  laminar  burning  velocities,  AL  is  the  flamelet  surface 
area  for  laminar  flow,  and  AT  is  the  ensemble  average  flamelet  area  in  turbulent  flow. 
Thus,  in  this  model  the  sole  effect  of  turbulence  is  to  generate  flamelet  surface  area. 

Damkohler’s  model  can  be  extended  to  higher  levels  of  turbulence  by  allowing  for 
perturbations  of  the  flamelet  structure  by  the  turbulence,  i.e.,  allowing  for  flamelet  stretch. 
Then 


ut =  <ul>  T*’  (2) 

L 

where  <uL>  is  the  mean,  stretched-flamelet  burning  velocity. 

Early  attempts  to  use  Damkohler's  model  suggestion  were  not  successful.  A  major 
reason  for  this  lack  of  success  appears  to  be  the  failure  of  area  ratio  estimates  to  account  for 
the  multiple  length  scale  wrinkling  of  flamelet  surfaces  by  turbulence  [5].  Rough  surfaces 
with  multiple,  self-similar  scales  of  wrinkling  can  be  characterized  with  the  mathematics  of 
fractals  [8].  In  our  research  we  have  used  this  characterization  to  estimate  the  area  ratio  in 
(1)  and  (2). 

Consider  an  isotropic,  homogeneous  fractal  surface  filling  a  volume  L3  in  space. 
One  aspect  of  the  surface's  fractal  character  is  that  if  one  measures  the  surface  area.  A,  in 

L3  with  a  measurement  scale  e  which  is  commensurate  with  the  scales  of  surface  wrinkling 

he  finds  for  varying  e  that 

A/L3~|e/L|'De2/L3.  (3) 

See  also  Fig.  1  and  [5, 8].  D  is  the  fractal  dimension  and  is  defined  by  (3).  Fractal 
character  is  reflected  in  the  power  law  dependency  of  A  versus  £. 

Mandelbrot  19]  has  suggested  that  constant  property  surfaces  in  isotropic, 
homogeneous  turbulence  are  fractal.  This  hypothesis  is  supported  by  experiment  [10,11] 


and  by  analysis  of  cloud  dispersion  [12].  These  surfaces  are  expected  to  exhibit  fractal 
character  for  a  limited  range  of  measurement  scales  corresponding  to  the  range  of  scales  of 
wrinkling.  The  experimental  results  of  [1 1]  indicate  that  the  limits  to  fractal  behavior  are 

associated  with  an  inner  cutoff,  which  is  the  Kolmogorov  (r|)  scale  of  turbulence,  and  an 
outer  cutoff,  which  is  the  integral  (/ )  scale.  This  finding  suggests  that  the  inertial  subrange 
turbulence  eddies  are  the  agents  of  self-similar  fractal  wrinkling.  A  simple  model  of  how 
an  eddy  could  wrinkle  a  surface  is  also  suggestive  of  this  conclusion.  An  eddy  of  length 
scale  /  n  and  velocity  scale  un  causes  a  sinusoidal  like  wrinkle  of  wavelength  /  n  and 

amplitude  u^At,  where  At  is  the  eddy  life  time.  Typical  scaling  gives  At  ~  /  „/  u„  and  the 
aspect  ratio  of  the  wrinkles  is  seen  to  be  one,  independent  of  the  wrinkle  amplitude  and 
wavelength,  ie.,  the  wrinkling  is  similar  at  all  scales. 

In  our  research  [5]  it  is  assumed  that  flamelet  surfaces  can  be  represented  as  fractal 
surfaces  and  the  assumption  is  used  to  estimate  the  area  ratio  in  (2).  It  is  assumed  that  there 

are  inner  (£t)  and  outer  (e^)  cutoffs  to  the  fractal  behavior  for  flamelets  and  that  the 
ensemble  average  of  measured  flamelet  area  will  vary  with  measurement  scale  as  shown  in 
Fig.  lb.  Further  it  is  argued  that  Ai,  the  area  at  the  inner  cutoff,  is  the  ensemble  average 
flamelet  area  in  L3,  while  A0,  the  area  at  the  outer  cutoff,  would  be  the  flamelet  area  in  the 
absence  of  turbulence.  Thus 


ut=  <UL>X’=<UL> 
0 


(4) 


One  sees  that  in  this  model  the  area  ratio  is  given  in  terms  of  three  parameters  -  two  length 
scales  and  the  fractal  dimension. 

To  complete  the  burning  velocity  model  expressions  for  these  length  scales,  a  value 
for  D  and  an  expression  for  <uL>  are  required.  In  [5]  it  is  argued  that 

e0  =  l,  and  (5) 

£j  =  T]  (l-(l-At‘1/4R/'3/4)  expf-A^R/-1'4  u’  /  <uL>  )).  (6) 


At  is  an  empirical  constant  obtained  from  data  on  noncombusting  turbulent  flow,  u'  is  the 
root  mean  square  of  velocity  fluctuations;  and  R/  is  the  turbulence  Reynolds  number  ( R /  = 

u7  /  v  ).  Also  presented  in  [5]  is  an  expression  for  <Ul>  based  on  an  analytical  expression 
for  flamelet  stretch  effects  developed  by  Law  [5].  For  brevity  this  expression  is  not 
repeated  here.  From  the  above  expressions  and  an  expression  for  the  ratio  1}  II  good  for 
grid  turbulence  one  finds 

ut  =  <uL>  {[  l-( l-At_1/4R/*3/4)  exp(-(A/R,)1/4  u'  /  <uL>  )]At1/4R/3'4}D-2.  (7) 

In  [5]  model  predictions  for  burning  velocity  are  compared  to  data  obtained  at 
Leeds  University  [13,  14].  A  value  for  D  of  2.37  is  used  based  on  the  results  of  [10-12]. 
The  comparison  shows  good  agreement  for  a  wide  range  of  R/  (1,000  -  40,000)  and  of 


u7u0  (  u7u0  £  40)  for  various  fuels. 


II.2.  Chemical  Closure  Model 

In  a  complementary  research  program  a  model  for  premixed  turbulent  flames  has 
been  developed  based  upon  conditional  statistics  and  a  second  order  transport  closure 
model  developed  by  Professor  Lumley  of  Cornell  and  his  students.  Model  equations  for 
conditional  second  order  statistics  are  proposed  and  solved.  Conditioning  is  on  being 
either  in  reactants  or  products  and  the  model  approach  assumes  combustion  to  occur  in 
flamelets.  A  paper  describing  this  model  will  appear  soon  [15]. 

We  have  extended  this  modeling  effort  as  part  of  our  ARO  supported  research  by 
proposing  a  chemical  closure  model  based  on  the  assumption  of  fractal  flamelets.  This 
work  was  started  at  Cambridge  University  (during  the  visit  of  F  C  G)  and  is  a  joint  effort 
with  Professor  KNC  Bray  of  Cambridge  and  Dr.  J-Y  Chen  of  Sandia  National 
Laboratories,  Livermore.  A  closure  model  has  been  proposed  for  the  mean  fuel 
consumption  rate,  and  it  has  been  evaluated  in  numerical  and  analytical  calculations. 

The  proposed  closure  model  expression  can  be  written  as 

<£of>  =  CRp0<AYfuL>f  (/  f/ri)D-2/  F-i  <c>  (1-<C>),  (8) 

with  f  given  by 

f  =  [l-(l-At'1/4R/'3/4)  exp(-At1/4R/-1/4u7<uL>)]. 

CR  is  a  model  constant;  p0  is  the  reactant  gas  density;  AYf  is  the  change  in  fuel  mass 
fraction  across  the  flamelet;  /  p  is  the  local  flame  brush  thickness;  and  c  is  the  reaction 
progress  variable  which  is  defined  as  the  ratio  of  the  fuel  mass  fraction  less  the  initial  fuel 

mass  fraction  and  AYf.  The  f  subscript  on  the  angle  bracket  denotes  a  conditional  average 
with  conditioning  on  the  flamelet. 

Details  of  the  model  development  and  the  analyses  are  presented  in  Appendix  A. 
Here  we  present  a  summary  of  the  findings.  The  model  is  used  in  simplified  analyses  of 
normal  and  oblique  flames. 

For  a  normal  flame,  we  assume  high  Reynolds  number,  one-dimensional  flow  in 

the  mean,  and  constant  AYf,  <uL>,  fl,f  and  / .  Integration  of  the  equation  for  <c>  across 
the  flame  brush  gives 

ut  /  <uL>  =  (/  f  /T))0*2 

(which  is  identical  to  our  burning  velocity  model  prediction)  provided  that 


which  is  reasonable.  Thus  the  closure  model  for  normal  flames  is  consistent  with  the 
burning  velocity  model. 

The  above  analysis  can  be  extended  to  the  oblique  flame  case  provided  several 
additional  assumption  are  made.  The  assumptions  to  be  made  include  constant  density, 
parabolic  flow  ( the  flame  brush  lies  almost  parallel  to  the  flow  direction  ),  and  flow 
similarity.  A  major  point  of  interest  in  this  analysis  is  to  see  if  a  turbulent  burning  velocity 
could  be  defined  for  the  oblique  flame  case.  A  definable  burning  velocity  is  obtained  with 

ut  /  <uL>  =  (/  f/r|)D*2  provided  (9)  is  valid  and  either  the  flames  brush  thickness  is  constant 
or  grows  linearly  with  distance  in  the  flow  direction. 

Numerical  calculations  using  the  transport  closure  of  [15]  were  performed  for 
oblique  flames  with  constant  density  in  decaying  grid  turbulence.  Calculated  unconditional 
normal  stresses  and  other  turbulence  properties  compared  well  with  the  grid  turbulence  data 
of  [16]  which  is  for  nonreacting  flow. 

Flame  results  were  obtained  for  a  range  of  conditions  with  a  single  value  for  CR;  CR 
=  4.0.  Calculated  burning  velocities  compare  favorably  with  values  obtained  from  the 
burning  velocity  model  as  would  be  expected  from  the  results  of  the  normal  and  oblique 
flames  analyses  discussed  above.  However  flame  brush  thickness  did  not  grow  linearly 
and  there  were  slight  but  noticeable  variances  between  the  burning  velocity  calculated  and 
that  obtained  from  the  ut  model.  These  differences  are  attributed  to  a  departure  from 
similarity  in  the  flows  obtained  by  numerical  calculations.  In  tum  the  departure  from 
similarity  is  attributed  to  turbulence  decay. 

As  part  of  the  numerical  calculations  we  obtained  results  using  the  new  chemical 
closure  model  with  a  two-equation,  gradient  transport  model  in  the  equation  for  <c>;  A  k-e 

model  was  used  with  k  and  e  obtained  from  the  second  order  model  results.  The 
differences  between  the  two  transport  models  are  striking.  The  gradient  model  under 
predicts  the  turbulent  fluxes,  and  consequently  the  flame  brush  thickness  is  also  under 
predicted.  On  the  other  hand  the  burning  velocities  predicted  using  the  new  chemical 
closure  model  and  either  of  the  two  transport  models  are  quite  similar. 

Overall  we  are  very  pleased  with  the  model  results.  The  representation  of  flamelet 
surface  geometry  by  fractals  appears  to  be  a  very  promising  approach  which  we  intend  to 
continue  to  pursue. 

II.3  Modeling  of  Jets  and  Jet  Flames 

The  flamelet  concept  may  also  be  applied  to  nonpremixed  flames  [  17].  In  this  case 
reaction  still  is  assumed  to  occur  in  thin  sheet  like  regions  which  are  centered  on  the 
stoichiometric  mixture  fraction  surface,  Z  =  Zs.  ( Z  is  the  mass  weighted  mixture  fraction. ) 
We  have  developed  modeling  ideas  applicable  to  jet  mixing  and  reacting  jets  based  on  the 
assumption  that  constant  mixture  fraction  surfaces  can  be  represented  as  fractal  surfaces.  A 
paper  applying  these  ideas  to  jet  mixing  has  been  submitted  for  journal  publication  [18], 
Appendix  B.  This  work  is  now  summarized. 


An  axisymmetric  jet  is  modeled,  and  the  jet  structure  is  viewed  as  a  set  of  constant 
concentration  surfaces.  Stationary  turbulence  and  uniform  density  are  assumed,  and  a 
balance  equation  for  the  ensemble  mean  flux  of  jet  fluid  across  a  Z  constant  surface  is 
developed  which  states  that  this  flux  must  equal  the  mean  flow  supplying  the  jet  Terms  in 
this  balance  equation  are  modeled  1)  by  assuming  that  the  geometry  of  constant 
concentration  surfaces  can  be  represented  by  fractal  surfaces  and  2)  by  estimating  the  mean 
jet  fluid  mass  flux  per  unit  area  of  a  constant  concentration  surface  in  terms  of  the 
Kolmogorov  velocity.  To  complete  the  analysis  similarity  in  the  jet  structure  is  assumed. 

Consider  a  slice  of  the  jet  defined  by  x  to  x  +  dx,  where  x  is  axial  distance.  The 
average  Z  surface  area  in  this  slice  is  modeled  as 

dAz  =  2jtr<2>(Al,/4R/3/4)D*2  dx  (10) 

r<z>  is  the  local  radius  of  the  <Z>  contour  with  <Z>  =  Z;  the  area  of  the  <Z>  constant 
surface  is  identified  as  A0.  The  average  flux  of  jet  fluid  per  unit  area  of  a  Z  constant 

surface  is  modeled  as  CdpZv  with  v  being  the  Kolmogorov  velocity,  p  the  density,  and  Cd 
a  model  constant.  The  average  flux  of  jet  fluid  across  a  Z  surface  between  x  and  x  +  dx  is 

obtained  by  multiplying  (10)  by  CdpZv.  The  resulting  expression  is  integrated  over  x 
using  a  standard  similarity  form  for  r<z>.  The  result  is 

L/d  =  C  Cd  R/  [3/(4-U)][D-2]  +  ll/2-3/(4-n)l  Zq.  (11) 

L  is  the  axial  distance  at  which  <Z>  on  the  jet  axis  falls  to  Zq,  and  d  is  the  initial  jet 

diameter.  C  is  an  empirical  constant  [18],  and  p  is  the  intermittency  exponent  --  0.25  £  p  <, 
0.5.  The  form  of  (1 1)  is  that  expected  from  similarity  provided  there  is  no  Reynolds 
number  dependency  which  in  turn  requires 


D  = 


,  Ox  +  2) 

2  +  - 5 - 


(12) 


Thus  the  jet  mixing  analysis  gives  an  expression  for  the  fractal  dimension.  The  expression 
is  identical  to  the  one  obtained  by  others  [12,  19]  from  an  analysis  of  cloud  dispersion  in 

homogeneous,  isotropic  turbulence.  For  the  expected  range  of  p  values,  (12)  gives  2.33  <, 
D  <,  2.42  which  is  consistent  with  experiment  and  the  value  used  in  our  burning  velocity 
model.  All  of  these  results  support  the  idea  of  a  universal  expression  for  D  valid  for  high 
Reynolds  number. 

Some  preliminary  work  has  been  performed  to  extend  our  fractal  modeling  ideas  to 
the  reacting  jet  case.  Variable  density  effects  have  been  added  in  an  approximation,  and  an 
expression  for  the  length  of  jet  flames  has  been  developed.  The  flame  length  expression 
with  D  satisfying  (12)  is 


L/d  =  C(l/Zs)^(p./pf) . 


(13) 


Zs  is  the  stoichiometric  jet  fluid  mixture  fraction,  while  pji and  Pf  are  the  initial  jet  fluid 
density  and  an  average  flame  gas  density,  respectively.  This  work  has  been  presented  and 
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a  report  is  available  [19]. 

As  in  the  case  of  premixed  flames  our  initial  application  of  fractal  concepts  to 
modeling  of  jet  mixing  and  combustion  has  proved  useful  and  encouraging.  Further  work 
with  fractals  on  nonpremixed  flames  is  planned  and  will  be  supported  by  ARO  in  a  new 
program  of  research  [20]. 


III.  EXPERIMENTAL  RESEARCH 


Our  experimental  research  may  be  divided  into  three  distinct  tasks:  measurement  of 
passage  times  in  turbulent  v-flames,  conditional  velocity  measurements  in  v-flames  and 
preliminary  studies  of  the  fractal  character  of  flamelets  in  low  Reynolds  number  turbulent 
premixed  flames.  The  results  of  each  of  these  efforts  are  reported  separately. 

III.l  Passage  time  measurements 

Bray,  Libby  and  Moss  (BLM)  have  proposed  a  model  for  premixed  turbulent 
flames  which  makes  use  of  the  probability  distribution  function  (pdf)  of  passage  times  -- 
the  time  duration  between  the  passage  of  a  flamelet  at  a  point  in  the  turbulent  flame  brush. 
In  their  work  BLM  distinguish  between  times  in  reactants  and  times  in  products.  In  their 
initial  work  [4]  the  pdf  was  assumed  to  take  the  form 

P(t)  =  -7-  exp(  -  —i-),  i  =  r  or  p.  (14) 

1  <tj>  <t.> 

r  and  p  refer  to  reactant  and  products  respectively,  while  <tj>  is  the  local  mean  passage 
time. 


This  form  is  not  satisfactory  at  small  tj  since  it  does  not  approach  zero  as  t* 
approaches  zero.  Recently  Bray  and  Libby  (BL)  have  proposed  another  form  for  P(tj) 

[21],  the  gamma  two  distribution. 

4  'i  2 

P(t.)  = - -exp(-— ■),  i  =  r,  p.  (15) 

<tj>  i 

This  distribution  approaches  zero  at  small  q  as  expected. 

We  have  made  Rayleigh  scattering  measurements  to  determine  P(ti)  and  test  the 
hypotheses  of  BLM  and  of  BL.  In  premixed  flames  Rayleigh  scattering  is  proportional  to 
the  gas  density,  and  where  flamelets  are  present,  the  temporal  signal  from  Rayleigh 
scattering  will  be,  apart  from  noise,  a  random  series  of  changes  from  high  to  low  and  back 
again.  The  changes  in  signal  level  correspond  to  a  flamelet  passage  event  and  therefore 
passage  time  data  can  be  obtained  from  Rayleigh  scattering  time  series  records  stored  on 
computer. 

Passage  time  measurements  were  performed  on  low  Reynolds  number  v-flames.  A 
schematic  of  the  burner  used  is  shown  in  Figure  1.  An  unconfined  v-flame  is  stabilized  in 
grid  turbulence  on  a  rod  (1.5mm  dia.)  mounted  across  the  exit  of  a  cylindrical  burner 
(50mm  dia.).  Turbulence  is  generated  by  woven-wire  screens  mounted  30mm  below  the 
rod.  Fuel  (commercial  grade  methane)  and  air  are  mixed  in  a  plenum  chamber  and  flow 
into  the  burner.  For  details  see  [2,  3]. 

Rayleigh  scattering  is  induced  by  a  1  watt  argon-ion  laser  operating  in  the  green. 
Scattered  photons  are  detected  by  a  photomultiplier  and  the  signal  is  stored  on  computer. 
The  spatial  resolution  of  the  measurements  is  approximately  0.1mm  and  the  temporal 

resolution  is  200|is.  Minimum  resolution  is  set  by  signal- to- noise  considerations  [3]. 

Scattering  measurements  were  performed  for  three  different  flame  conditions  and  at 
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several  different  spatial  locations  in  each  flame.  The  initial  results  for  the  three  flames, 
which  are  reported  in  [3],  are  in  rough  agreement  with  the  hypothesis  of  BLM  at  least  for 
large  tr  There  are  discrepancies  in  the  comparisons  at  small  tj  where  the  gamma  two 
distribution  suggested  by  BL  is  expected  to  give  a  better  result  Comparison  of 
experimental  data  with  a  gamma  two  distribution  does  show  better  agreement  for  the  lower 
Reynolds  number  flames  at  smaller  q.  But  comparisons  at  larger  tj  and  for  the  higher 
Reynolds  number  flame  are  not  as  good  as  given  by  (14).  Sample  pdfs  are  shown  in  Fig. 

2.  On  the  basis  of  these  data  one  concludes  that  there  is  room  for  further  improvement  in 
the  representation  of  P(t,). 

III.2  Conditional  Velocity  Measurements 

For  low  Reynolds  number  turbulence,  as  noted  above,  the  flamelet  structure  is 
expected  to  approach  that  of  an  unperturbed  laminar  flame.  Assuming  an  unperturbed, 
quasi-steady  structure  there  is  a  jump  in  fluid  velocity  normal  to  the  flamelet  surface  across 
this  surface  which  is  equal  to  the  laminar  burning  velocity  times  the  density  ratio  across  the 
flamelet  minus  one.  For  small  u'  relative  to  u0  this  acceleration  could  be  important  to 
turbulence  generation.  We  undertook  conditional  velocity  measurements  in  order  to  detect 
flow  acceleration  across  flamelets.  These  measurements  entailed  the  simultaneous 
measurement  of  flamelet  position  and  gas  velocity.  Laser  induced  scattering  techniques 
were  employed  in  these  measurements.  Some  of  the  results  have  been  read  at  a  ASME- 
JSME  meeting  and  are  available  in  report  form  [7]  and  Appendix  C. 

These  measurements  were  made  on  v-flames.  Gas  velocity  is  measured  with  a 
commercial  LV  system  composed  of  an  argon-ion  laser,  counter-type  signal  processor,  and 
associated  optical  components.  The  spatial  resolution  of  the  velocity  measurements  is  a 
few  tenths  of  a  millimeter. 

Flamelet  position  alone  a  line  in  space  is  measured  with  a  light  scattering  technique. 
The  reactant  flow  is  seeded  with  an  oil  mist  which  evaporates  in  the  flamelet  and  thus 
marks  reactants  but  not  products.  Scattering  from  a  He-Ne  laser  is  detected  with  a 
photomultiplier  tube  which  sees  scattering  from  a  one  centimeter  portion  of  the  laser  beam. 
Thus  the  detected  signal  strength  depends  on  the  amount  of  reactant  gas  along  the  portion 
of  the  laser  beam  viewed,  i.  e.,  the  position  of  the  flamelet  along  the  laser  beam.  The  He- 
Ne  laser  beam  is  alined  in  the  horizontal;  see  Fig.  1. 

For  measurements,  the  LV  measurement  volume  is  located  at  different  points  along 
the  He-Ne  laser  beam.  Position  and  velocity  data  are  recorded  simultaneously  on  a 
computer.  The  oil  mist  used  for  the  position  measurements  is  excellent  for  LV 
measurements.  However,  since  the  oil  evaporates  in  the  flamelet,  a  second  refractory 
particle  seed  must  be  added  to  the  flow  for  velocity  measurements  in  the  products.  The 
refractory  seed  levels  are  low  enough  so  that  the  position  measurement  is  not  affected  by 
the  addition  of  refractory  scattering  particles.  Vertical  and  horizontal  (perpendicular  to  the 
vertical  plane  containing  the  stabilizer  rod)  velocities  are  measured. 

The  position  and  velocity  data  are  analyzed  to  give  the  pdf  of  flamelet  position  along 
the  He-Ne  laser  beam  and  conditioned  velocities.  Conditioning  is  by  being  either  in 
products  or  reactants  and  by  distance  from  the  flamelet.  We  expect  that  conditioning  by 
distance  will  provide  evidence  regarding  flow  acceleration  across  flamelets. 

Measurements  in  reactants  (without  refractory  seed)  have  been  performed  and 
results  reported  [7].  Measurements  with  refractory  seed  are  in  progress.  They  are  the 
thesis  research  of  Paul  Miles;  he  is  expected  to  complete  his  thesis  in  August  1987.  Results 
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to  date  of  note  are  summarized  in  the  following. 

The  pdf  of  flamelet  position  is  nearly  Guassian  and  the  its  standard  deviation  is  a 
useful  measure  of  the  flame  brush  thickness.  The  standard  deviation  grows  approximately 
as  a  linear  function  of  downstream  distance.  The  derivative  of  the  flamelet  position  signal 
is  determined  and  viewed  as  a  measure  of  the  flamelet  velocity  [7].  Surprisingly  high 
values  of  the  derivative  are  observed.  Power  spectra  of  the  derivative  signal  are  compared 
to  spectra  of  the  horizontal  velocity  (the  V  component).  The  two  spectra  are  nearly 
coincidence  in  the  frequency  range  from  approximately  500  to  about  1500  Hz.  Mean 
velocities  measured  in  reactants  are  consistent  with  those  measured  by  others  in  a  similar  v- 
flame.  Velocity  statistics  conditioned  by  distance  from  the  flamelet  in  reactants  do  not 
contain  surprises.  In  particular  rms  levels  do  not  vary  significantly  with  distance  —  see 
Fig.  10,  Appendix  C. 

Measurements  with  refractory  seed  are  in  progress.  The  persistence  due  to  slow 
evaporation  of  large  oil  drops  through  the  flamelet  have  caused  problems  in  these 
measurements.  Refractory  seed  concentrations  are  orders  of  magnitude  below  that  of  the 
oil  mist.  As  a  consequence  a  small  percentage  of  large  oil  drops  which  do  not  influence 
the  position  measurement  or  the  reactant  velocity  measurement,  appear  to  bias  the  product 
gas  velocity  measurements.  We  are  currently  working  to  resolve  this  problem. 

III.3  Evaluation  of  the  fractal  character  of  flamelets 

While  it  is  not  obvious  how  to  assess  surface  fractal  character  directly,  inferences 
can  be  draw  from  the  application  of  laser  tomography  and  point,  scalar  measurements;  see 
for  example  [11].  Preliminary  experiments  on  v-flames  have  been  performed  using 
tomography  and  point  scattering  measurements  [24, 25]. 

The  intersection  of  a  flamelet  surface  with  a  plane  can  be  visualized  with  laser 
tomography  [23].  For  these  experiments  the  reactant  flow  is  seeded  with  an  oil  mist  as  for 
the  velocity  measurements.  The  output  of  a  high  energy,  pulsed  Nd:Yag  laser  is  frequency 
doubled  into  the  green  (300  mj  per  pulse  in  the  green),  and  lenses  used  to  form  the  laser 
beam  into  a  sheet  of  light  approximately  0.5mm  in  thickness.  The  flame  brush  is 
illuminated  by  the  laser,  and  the  reactant  gas  in  the  plane  of  illumination  is  made  visible  by 
scattering  of  light  from  the  oil  mist.  Photographs  of  the  flame  brush  are  taken  with  a  35mm 
camera  synchronized  with  the  laser  pulse.  Because  the  oil  evaporates  in  the  flamelet  there  is 
no  scattering  from  the  products,  and  consequently  the  curve  defined  by  the  boundary 
between  regions  of  scattering  and  of  no  scattering  observed  in  the  photographs  -called 
tomograms  -  is  the  intersection  of  the  plane  defined  by  the  laser  and  the  flamelet  surface. 

Curves  formed  by  the  intersection  of  a  plane  with  a  fractal  surface  are  fractal,  and 
for  isotropic  fractal  surfaces  the  fractal  dimension  of  the  curve  is  one  less  than  that  of  the 
surface  from  which  it  is  formed.  We  assume  that  the  curve  fractal  dimensions  which  we 
measure  are  also  one  less  than  that  of  the  fractal  flamelet  For  curves,  fractal  dimension  can 
be  defined  in  terms  of  measured  length,  L ,  in  a  manner  analogous  to  that  for  surfaces. 

Thus 


L  /  L2  ~  |  e  /  L  |'De  /  L2.  (16) 

Here  D  is  the  fractal  dimension  of  the  curve;  L  is  the  scale  of  the  surface  area  over  which  L 
is  measured;  and  e  is  the  measurement  scale. 

In  our  studies  tomograms  for  several  different  flame  conditions  —  see  Table  1  —  are 


digitized  and  their  fractal  dimension  determined  by  application  of  (16).  Plots  of  L  versus  e 
were  made  and  slopes  taken.  A  sample  plot  is  given  in  Fig.  3.  Since  there  are  limits  to  the 
length  scales  of  wrinkling  we  expect  a  curve  similar  in  shape  to  that  in  Fig.  lb,  two 
horizontal  lines  connected  by  a  straight  line  with  slope  1  -  D.  Furthermore  for  a  surface 
fractal  dimension  of  2.37  we  expect  a  curve  dimension  of  1.37.  In  Fig. 3  we  see  that  the 

shape  is  in  general  what  we  expect.  However  the  maximum  e  is  not  large  enough  to 
unambiguously  show  behavior  above  the  outer  cutoff.  Also  the  fractal  dimension  taken 
from  the  curve  is  well  below  the  expected  value.  Low  fractal  dimensions  are  observed  for 
all  the  flame  conditions  studied,  Table  1.  From  measured  length  curves,  inner  and  outer 
cutoffs  are  also  estimated  and  found  to  be  several  times  larger  than  the  Kolmogorov  and 
integral  scales  respectively. 

We  attribute  the  low  fractal  dimension  obtained  in  these  measurements  to  two 
factors,  low  turbulence  Reynolds  number  and  low  u'/uq.  Fractal  wrinkling  is  most  likely 
caused  by  the  inertial  subrange  turbulent  eddies.  For  low  Reynolds  number  turbulence 
these  eddies  are  not  fully  developed  which  may  influence  the  character  of  surface 
wrinkling.  In  section  II.  1  a  simple  model  for  surface  wrinkling  by  turbulent  eddies  was 

presented.  For  this  model,  the  amplitude  of  wrinkling  depends  on  a  time  scale.  At,  which 
is  the  time  a  surface  is  exposed  to  an  eddy  and  was  taken  as  the  eddy  life  time.  Since  a 
flamelet  propagates  relative  to  the  reactant  mixture,  exposure  time  may  be  less  than  an  eddy 
life  time.  An  alternate  exposure  time  is  given  by  ln  /  u0.  For  u0  >  u„  this  latter  time  is  the 
exposure  time  and  the  amplitude  of  corrugation  by  our  model  is  /  „  (Uj/uq)  Thus  the 
condition  un  =  u0  marks  a  transition  in  behavior.  The  /  n  for  un  =  u0  is  defined  as  the 
Gibson  scale  -- 1  q~  and  by  standard  scaling  arguments  /  G  =  /  (uq  /  u’  )3.  For  our 
conditions  /  G  >  / .  By  our  simple,  wrinkling  model  all  wrinkle  amplitudes  scale  as  /  n 
(un/uQ)  rather  than  l  n  and  it  is  reasonable  to  expect  a  variation  in  fractal  behavior.  By  the 
above  argument  one  expects  D  would  approach  1.37  for  large  R/  and  large  u'/uq. 

The  set  of  points  formed  by  the  intersection  of  a  line  with  a  fractal  surface  also 
exhibits  fractal  character,  and  for  an  isotropic  fractal  surface  the  fractal  dimension  of  the  set 
of  points  is  two  less  than  the  fractal  dimension  of  the  surface.  A  fractal  set  of  points  is 
known  as  a  fractal  dust  [8],  For  points  along  a  line  the  fractal  dimension  is  determined  by 
the  variation  in  the  measure  of  the  set  with  measurement  scale.  Divide  the  line  into 

segments  of  length  e.  The  number  of  segments  in  which  there  is  one  or  more  members  of 
the  set  is  N.  Then  for  a  fractal  set 


N  "  £*d.  (17) 

D  is  again  the  fractal  dimension,  and  its  value  falls  between  0  and  1. 

Time-series,  point-scattering  measurements  were  performed  to  study  the  fractal 
character  of  the  distribution  of  flamelet  crossing  events  in  time.  As  before  the  reactants  are 
seeded  with  an  oil  mist  and  scattering  is  excited  by  a  He-Ne  laser.  Scattered  photons  from 
a  small  volume  are  detected  and  used  to  drive  a  logic  circuit  which  generates  a  short 
duration,  logic  pulse  at  each  flamelet  crossing  event.  These  pulse  control  a  computer  clock 
which  counts  the  time  between  crossing  events,  and  these  times  are  stored  in  computer. 
The  spatial  resolution  of  these  measurements  is  a  few  tenths  of  a  millimeter  and  the 

temporal  resolution  is  50  (isec. 
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From  recorded  data  the  measure  of  a  set  of  crossing  events  along  the  time  axis  can 
be  determined  as  a  function  of  e.  Plots  of  N  versus  e  are  then  generated;  see  Fig.  4. 

Fractal  behavior  is  not  obvious  in  the  data  shown  in  Fig.  4.  These  data  may  be  represented 
by  several  straight  line  segments  and  a  physical  interpretation  can  be  given  to  each  of  the 

segments.  The  horizontal  line  to  the  left  is  for  the  case  where  e  is  below  the  inner  cutoff 
and  all  members  of  the  set  are  resolved.  The  right  most  line  has  a  slope  of  -1  which 

indicates  e  values  above  an  outer  cutoff.  In  each  segment  there  are  set  members  for  all  £  in 
this  range,  and  none  of  the  time  intervals  between  crossing  events  are  resolved  for  these 

large  £'s.  In  between  these  two  regions  of  limiting  behavior  are  two  regions  of  fractal 
behavior;  one  having  a  fractal  dimension  of  approximately  0. 1  and  the  other  a  dimension  of 
approximately  0.5.  The  tomographic  results  imply  a  fractal  dimension  for  the  dust  of 
approximately  0.1,  and  therefore  we  associate  this  segment  of  the  curve  with  the  fractal 
wrinkling  of  the  flamelet  surface. 

To  understand  the  portion  of  the  curve  where  D  ~  0.5  consider  a  coin  toss  game  in 
which  the  winner  of  each  toss  receives  the  coin.  The  coins  are  identical  and  fair.  A  plot  of 
cumulative  winnings  for  one  of  the  players  versus  trial  number  is  a  random  walk,  while  the 
distribution  of  zero  winnings  versus  trial  number  is  a  fractal  set  with  dimension  1/2  [8]. 
This  observation  suggests  that  the  D  -  0.5  segment  may  be  the  result  of  large  scale  flapping 
of  the  flamelet.  In  other  words,  a  smooth,  flapping  flamelet  would  give  fractal  behavior 
with  D  -  0.5,  while  the  fractal  wrinkling  of  the  flamelet  contributes  a  region  in  our  N 

versus  e  plot  with  D  ~  0. 1.  With  regard  to  the  proposed  burning  velocity  and  chemical 
closure  models  the  fractal  character  and  fractal  dimension  of  the  surface  wrinkling  are  the 
critical  features  to  be  tested  by  the  experiments.  Therefore  we  conclude  tentatively  that  the 
experiments  support  the  concept  of  fractal  wrinkling  and  for  the  conditions  studied  a  fractal 
dimension  of  approximately  2.1.  We  speculate,  based  on  experiments  [10, 1 1]  and 
analyses  for  nonreacting  flows  [12,  18,  19]  that  as  u’/u0  and  R/  increase  D  will  approach 
the  limiting  value  given  by  (12). 

Our  investigations  of  fractal  behavior  are  a  preliminary,  first  step,  and  our 
interpretation  of  results  is  also  preliminary.  More  work  is  required,  especially  at  higher 
Reynolds  number. 


IV.  TURBULENCE  GENERATION 

For  the  flames  we  have  studied  reaction  occurs  in  thin  sheets  "separating  low  and 
high  temperature  regions  of  uniform  density.  From  a  mechanistic  point  of  view  turbulence 
generation  associated  with  density  change  in  these  flames  must  be  associated  with  the 
flamelet.  In  unconditioned  moment  closure  models  density  fluctuation  effects  appear  in 
many  terms  which  require  modeling.  However  these  terms  do  not  explicitly  account  for  the 
fact  that  density  fluctuations  effects  are  found  only  in  the  flamelet.  We  believe  that  this  is  a 
serious  short  coming  of  these  models  because  a  more  detailed,  mechanistic  picture  of 
density  change  effects  when  obtained  from  experiment  (  or  other  sources )  cannot  be 
incorporated  explicitly  in  the  model.  On  the  otherhand  such  mechanistic  information  can  be 
incorporated  in  other  model  approaches  such  as  the  conditional  approach  developed  at 
Cornell  [15]. 

Similar  comments  apply  to  the  modeling  of  chemical  species  transport  effects  in 
premixed  flames  dominated  by  flamelets.  Temperature  and  composition  are  nearly  uniform 


in  reactants  and  products.  Molecular  transport  of  heat  and  species  occurs  in  the  flamelet. 
Unconditioned  moment  closures  cannot  reflect  these  facts  directly  while  other  model 
approaches  can. 

The  conditioned  velocity  measurements  were  undertaken  primarily  to  investigate 
flow  acceleration  across  the  flamelet.  Unfortunately  the  persistence  of  oil  droplets  across 
the  flamelet  hinders  a  quantitative  evaluation  of  such  acceleration.  We  continue  to  work  on 
this  problem. 


V.  SUMMARY 

Several  tasks  of  research  on  turbulent  flames  have  been  pursued.  The  focus  has 
been  on  low  to  moderate  Reynolds  number  premixed  flames.  A  fractal  representation  of 
surfaces  in  turbulent  flows  is  found  to  be  a  very  useful  tool  for  modeling  purposes. 
Experiments  support  the  hypothesis  that  flamelet  surfaces  can  be  represented  by  fractal 
surfaces,  but  more  experiments  are  need  before  the  evidence  can  be  considered  conclusive. 
Passage  time  pdfs  are  also  useful  in  combustion  modeling.  Such  distributions  were 
measured  and  found  to  be  reasonably  well  represented  by  one  of  two  proposed  functional 
forms  for  the  pdfs.  If  a  very  exact  representation  is  required  neither  function  seems  ideal. 
Flame  position  and  conditioned  velocity  measurements  were  performed.  The  pdf  of 
flamelet  position  is  nearly  guassian  and  the  standard  deviation  of  the  distribution  is  a 
measure  of  flame  brush  thickness.  The  velocity  measurements  reveal  the  flamelet  to  be  a 
source  of  flow  acceleration.  But  experimental  problems  have  hindered  the  gathering  of 
quantitative  data  on  the  magnitude  of  this  acceleration. 
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Table  1.  Tomographic  Fractal  Analysis  Results 
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8  mesh  screen  turbulence  generator 


$ 


0.8 

1.0 

0.8 


mean  axial 

number  of  curves 

mean  fractal 

standard 

velocity  (m/s) 

JU 

analyzed 

dimension 

deviation 

4.1 

25 

18 

1.09 

0.03 

4.2 

25 

16 

1.11 

0.03 

3.8 

23.5 

16 

1.11 

0.04 

1  7 

FIGURE  CAPTIONS 

1 .  Schematic  of  burner  showing  location  of  wire  mesh  turbulence  generator,  stabilizer 
rod,  and  laser  beam  for  flamelet  position  measurement.  The  inner  jet  diameter  is  5 
cm,  and  the  outer,  annular  jet  diameter  is  7.6  cm.  The  fuel  is  methane. 

2.  Distribution  of  passage  times  in  the  10  mesh,  <U>  =  3.8  m/s  flame.  Squares 
denote  times  in  products  and  pluses  denote  times  in  reactants  at  a  point  5.5  cm 
downstream  from  the  wire  mesh.  Triangles  denote  times  in  products  and  X's  times 
in  reactants  at  7.0  cm.  <c>  =  0.46.  The  straight  line  indicates  the  exponential  pdf 
torm  while  the  curved  line  is  the  gamma  2  distribution.  See  Ref.  3  for  details. 

3 .  Typical  length  plot  for  fractal  analysis  of  tomograms  [24]. 

4 .  Typical  plot  for  fractal  dust  analysis  of  the  time  series  measurements  of  flamelet 
crossing  events  [25]. 
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Abstract 


A  chemical  closure  model  for  premixed  turbulent  flames  is  proposed  and  tested  by 
analysis  and  numerical  computation  for  flames  with  vanishingly  small  density  change.  The 
model  is  based  on  the  assumption  that  reaction  zones  can  be  modeled  as  thin  sheets  — 
flamelets  --  and  that  the  geometry  of  these  sheets  can  be  represented  by  fractal  surfaces. 
The  model  expression  for  mean  fuel  consumption  rate  is 

<cof>  =  CRp0  <AYfUL>f  (/  f/  ti)0-2/  p-1  <c>(l-<c>), 

with  f  given  by 

f  =  [1  -  (1  -  At‘1/4R/‘3/4)exp(-At1/4R/'1/4u'/<uL>f)], 
and 

/  /  "H  =  Atl/4Rf 3/4, 

where  D  is  the  fractal  dimension  of  the  flamelet  surface  and  is  the  new  parameter  introduced 
by  the  fractal  geometry  assumption. 

This  model  is  tested  in  simplified  analyses  of  normal  and  oblique  flames  with  good 
results.  The  oblique  flame  analysis  provides  new  insight  into  the  definition  of  the  turbulent 
burning  velocity.  Numerical  computations  are  performed  with  a  conditioned  second  order 
closure  scheme,  and  the  chemical  closure  model  performance  is  found  to  be  good. 
Computed  results  with  a  gradient  transport  model  for  species  diffusion  show  that  turbulent 
fluxes  are  significantly  under  predicted  in  comparison  with  the  second  order  closure 
results. 


I.  Introduction 


In  [1]  a  model  for  turbulent  burning  velocity  is  proposed  based  on  the  assumption 
that  chemical  reaction  and  heat  release  rates  can  be  modeled  as  occurring  in  thin  sheets  — 
flamelets  --  and  that  the  geometry  of  these  sheets  can  be  represented  as  fractal.  It  is  argued 
that  flamelet  surfaces  are  extremely  rough  with  multiple  scales  of  wrinkling  and  that 
mathematical  relationships  true  for  fractal  surfaces  can  be  used  to  advantage  in 
characterizing  flamelet  surface  area. 

Fractal  surfaces  are  rough  and  a  significant  feature  of  this  roughness  is  the  regular 
variation  of  measured  surface  area  with  measurement  scale  [1,2,3]-  For  an  isotropic  fractal 
surface  in  a  volume  L3, 

A  /  L3  ~  e2  D  /  L3.  (1) 

A  is  the  measured  area,  and  e  is  the  measurement  scale.  D  is  the  fractal  dimension  of  the 
surface  and  is  bounded  such  that  2  <D  <  3.  Fractal  surfaces  are  extremely  wrinkled,  and 
according  to  Mandelbrot  [3]  they  appear  to  be  space  filling.  The  larger  D  the  more  space 
filling  a  fractal  surface  will  appear. 

For  flamelets  as  for  constant  property  surfaces  in  turbulent  flow,  there  is  a  limit  to 
the  range  of  length  scales  over  which  fractal  behavior  is  observed.  The  variation  of  A  with 
e  will  not  follow  (1)  for  e  outside  the  fractal  range,  and  for  surfaces  with  a  limited  range  of 
fractal  behavior  Gouldin  [1]  hypothesizes  a  simple  behavior  for  A  versus  e,  which  is 

depicted  in  Fig.  1.  At  this  point  two  new  parameters  enter  the  fractal  description,  the  inner 
and  outer  cutoffs  --  e*  and  e0  [3].  Flamelet  surface  area  is  finite  and  thus  A  must  approach 
a  limit ,  Aj,  as  £  -»  0.  For  large  e  and  anchored  flames  such  as  a  bunsen  flame,  we  expect 
that  on  large  measurement  scales  the  surface  will  appear  smooth  and  A  — >  A0  as  e  — » °°. 


This  limiting  behavior  at  both  large  and  small  e  is  represented  by  horizontal  line  segments 
with  £,  and  £0  defined  such  that  A0/  L3  ~  £02'D  /  L3  and  A;  /  L3  ~  Ei2‘D  /  L3. 

Mean  flamelet  surface  area  is  estimated  in  [1]  using  the  results  displayed  in  Fig.  1. 
Here  we  use  similar  reasoning  to  that  in  [1]  to  obtain  a  model  for  the  local  mean  rate  of  fuel 
consumption.  Then  normal  and  oblique  flames  are  analyzed  with  the  proposed  rate  model 
under  the  assumption  that  density  change  is  negligible,  and  a  burning  velocity  expression  is 
obtained  which  is  identical  to  that  of  [1],  A  numerical  analysis  of  anchored  flames  in  grid 
turbulence  is  performed  with  the  new  closure  model  and  with  a  conditioned,  second  order 
transport  model  developed  by  Chen,  Lumley  and  Gouldin  [4].  The  calculation  results  are 
realistic  and  consistent  with  those  of  [1].  Finally,  flame  calculation  results  with  a  gradient 
transport  model  for  turbulent  species  flux  are  compared  to  results  from  the  second  order 
model.  The  gradient  model  predictions  are  not  in  agreement  with  those  of  the  second  order 


II.  Chemical  Closure  Model 


The  proposed  rate  expression  for  fuel  consumption  is  of  the  form 
<cof>  =  CRp0  <AYfuL>f  (ej/e0)2-D/  c~l  <c>( l-<c >)/  c//  F ,  (2) 

where  CR  is  a  model  "constant",  Ul  is  the  local  laminar  burning  velocity  of  the  flamelet  and 
AYf  is  the  change  in  fuel  mass  fraction  across  the  flamelet.  The  notation  <  >f  represents  an 
average  over  the  flamelet  surface,  /  p  is  the  local,  turbulent  flame  brush  thickness,  while  /  c 
will  be  defined  below.  Finally,  p0  is  the  cold  gas  density,  and  <c>  is  the  mean  degree  of 
reaction  where  c  =  (Yf  -  Yf°)/AYf  and  Yf°  is  the  initial  fuel  mass  fraction.  This  expression 

is  justified  using  the  following  arguments. 

Consider  a  cubic  volume  centered  in  space  on  a  point  x  with  side  /  c  which  is  small 

enough  that  the  ensemble  mean  conditions  inside  the  volume  can  be  assumed  to  be  uniform 
with  a  mean  degree  of  reaction  <c>.  Assume  that  the  flamelet  geometry  can  be  represented 
by  a  fractal  surface  of  fractal  dimension  D,  inner  cutoff  Ej,  and  outer  cutoff  e0.  Here,  e0  is 

representative  of  the  larger  scales  of  wrinkling  of  the  flamelet  surface  and  for  the  flamelet 
surface  to  be  distributed  on  average  in  a  uniform  manner  over  /  c3  we  require  that  /  c  £  £0. 
Let  Pc  be  the  probability  that  a  flamelet  is  present  in  the  volume  /  c3. 

According  to  Mandelbrot  [1,3]  a  constant  property  surface  (  e.g.,  a  constant 
temperature  surface)  in  homogenous,  isotropic  turbulence  is  fractal.  In  a  cube  of  volume 
L3,  the  mean  surface  area  per  unit  volume  of  a  fractal  surface  measured  at  a  length  scale  e 

is 

A/L3  =  C|e/L|'De2/L3. 

C  is  a  constant  termed  by  Mandelbrot  the  lacunarity  [2].  As  noted,  Gouldin  [1]  has  argued 
that  if  there  is  a  characteristic  minimum  length  scale  of  wrinkling,  which  is  defined  as  the 
inner  cutoff,  the  actual  mean  surface  area  per  unit  volume  is  finite  and  can  be  written  a 
A/L3  =  C  |£j/L  |'d£j2/L3. 
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Mandelbrot  considered  homogeneous,  isotropic  turbulence.  Then  L3  is  a  volume 
embedded  in  turbulent  fluid  and  the  surface  is  dispersed  throughout  this  volume  in  a 
statistically  uniform  manner.  Flamelets  are  not  expected  to  fill  space  in  the  same  manner  as 
Mandelbrot's  constant  property  surfaces.  To  adapt  Mandelbrot's  fractal  surface  concepts 
to  the  present  problem  it  is  assumed  that,  when  there  is  flamelet  present  in  /  c3,  it  fills  /  c3 

in  the  same  manner  as  a  constant  property  surface  in  isotropic,  homogeneous  turbulence 
fills  L3.  Thus  the  average  flamelet  surface  area  in  /  c3  can  be  expressed  as 

Ac  =  C|ei//c|-D£i2Pc,  (3) 

where  Pc,  as  noted,  is  the  probability  of  finding  a  flamelet  in  the  volume  /  c3. 

For  modeling  we  propose  that  the  mean  rate  of  fuel  consumption  can  be  written  as  a 
product  of  Ac  and  the  mean  rate  of  fuel  consumption  per  unit  flamelet  area  --p0<AYfUj_>f. 

Flamelet  surface  wrinkling  and  other  effects  of  turbulence  will  perturb  the  flamelet  structure 
and  therefore  alter  fuel  consumption  rates.  p0AYfuL  represents  the  local  instantaneous 

consumption  rate  in  terms  of  the  reactant  density,  the  change  in  fuel  mass  fraction  across 
the  flamelet  and  the  laminar  burning  velocity.  This  expression  may  be  interpreted  as  a 
definition  for  the  laminar  burning  velocity.  It  is  implicit  in  our  modeling  approach  that 
good  estimates  of  p0<AYfUL>f  can  be  obtained  and  our  focus  is  on  the  modeling  of  Ac.  A 
fuller  discussion  of  possible  modeling  for  p0<AYfuL>f  is  given  in  [1]. 

With  our  expression  for  mean  fuel  consumption  per  unit  area,  the  mean 
consumption  rate  in  /  c3  is  just  p0<AYfUL>fAc  and 
<cuf>  ~  p0<AYfUL>f  et2  D/  c^Po 
with  Pc,  ex  and  /  c  yet  to  be  determined. 

There  is  experimental  evidence  for  v-flames  [5]  that  the  probability  density  function 
describing  flamelet  location  on  a  line  which  is  normal  to  the  flame  brush  is  Gaussian.  If 
flamelet  motion  is  the  result  of  many  random  impulses  the  central  limit  theorem  also 
suggests  that  this  probability  will  be  Gaussian.  For  thin  flamelets  a  reasonable 
approximation  [5]  to  the  value  of  the  Gaussian  function  at  a  point  x  is  given  by 
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Cn<c(x)>(l-<c(x)>)  where  Cn  is  a  normalization  constant.  For  clarity  the  functional 

dependency  of  <c>  on  position,  x,  is  shown  explicitly.  The  probability  of  Finding  the 
flamelet  along  a  segment  of  the  normal  which  is  /c  is  just  Cn  <c>(l-<c>)  l  c.  The 
integrated  probability  over  the  flame  brush  thickness,  l  p,  should  be  one  and  therefore  Cn  - 
/  p.  For  Pc  it  follows  that 

Pc  ~  <o(l-<c>)/c//F  .  (4) 

The  inner  cutoff  is  a  measure  of  the  smallest  scales  of  wrinkling.  For  unit  Schmidt 
number  we  expect  that  the  smallest  possible  scale  of  wrinkling  due  to  turbulent  velocity 
fluctuations  is  related  to  the  Kolmogorov  microscale.  Gouldin  [1]  has  hypothesized  that 
for  flamelets  the  inner  cutoff  is  equal  to  or  greater  than  the  Kolmogorov  scale  and  that  it 
depends  on  the  ratio  u'/<UL>f.  He  suggests  an  expression  for  ejl]  - 

£i  =  rj/f  =  r\/[l  -  (1  -  At’1/4R/'3/4)exp(-Atl/4R/’1/4u7<uL>)],  (5) 
where  tj  is  the  Kolmogorov  microscale,  R/  the  turbulent  Reynolds  number  and  u'  is  the 
root  mean  square  of  velocity  fluctuations.  T1  represents  the  minimum  possible  scale  of 
surface  wrinkling,  while  f  accounts  for  smoothing  at  the  smallest  scales  of  wrinkling  due  to 
flamelet  propagation  relative  to  the  reactant  flow.  While  this  relationship  for  f  is  speculative 
and  future  research  may  well  lead  to  an  alternate  expression  for  e,,  the  above  form  is 

adopted  here. 

The  arguments  in  [1]  suggest  that  /  c  should  be  a  multiple  of  / ,  the  turbulence 

integral  scale,  if  not  equal  to  l .  Implicit  in  the  arguments  leading  to  (3)  is  the  assumption 
that  when  present  the  wrinkled  flamelet  surface  is  assumed  to  fill  /c3  as  Mandelbrot's 
fractal,  constant  property  surface  fills  L3.  In  order  for  this  to  be  true  we  expect  /«.»/, 
where  /  is  taken  as  the  outer  cutoff  (/  =  ej.  Our  argument  here  is  based  on  the  assumption 
that,  while  the  flamelet  can  extend  over  large  distances  in  two  directions,  its  extent  in  the 
third  direction  is  limited.  At  any  instant  flamelet  sheets  can  be  bounded  by  two  smooth 
surfaces  with  a  local  separation  distance  which  is  proportional  to  the  outer  cutoff  which  in 
turn  is  assumed  to  be  /.  If  flamelet  wrinkling  is  such  that  flamelets  double  back  on 


themselves  many  times  the  separation  distance  of  these  two  bounding  surfaces  would  be 
larger  and  /c  could  be  larger  than  the  outer  cutoff.  We  assume  that  such  doubling  back  on 

large  scales  is  not  important.  Doubling  back  on  scales  up  to  the  order  of  the  outer  cutoff  is 
expected  and  accounted  for  by  the  fractal  description  of  flamelets 

At  this  point  a  digression  on  the  length  scales  introduced  in  this  work  may  be 
helpful.  r\  and  /  are  traditional  turbulence  length  scales  and  are  widely  used  for 
characterizing  turbulence.  £j  and  E0  are  scales  introduced  with  the  fractal  description  of 

rough  surfaces,  and  they  are  related  by  hypothesis  (which  is  supported  by  experiment)  to 
the  turbulence  length  scales,  r\  and  / .  /  c  is  introduced  for  purposes  of  exposition. 

A  major  difficulty  in  adapting  Mandelbrot's  results  [3]  for  surfaces  in  turbulent 
flow  to  the  present  problem  arises  becacuse  Mandelbrot  assumes  isotropic  surface 
distributions  while  our  flamelet  distributions  are  highly  anisotropic.  To  overcome  this 
problem  we  focus  attention  on  a  small  volume,  /c3.  When  the  flamelet  is  in  this  volume  its 

distribution  in  the  volume  is  nearly  isotropic,  and  an  ensemble  of  such  occurrences  should 
have  an  even  more  isotropic  distribution.  Clearly,  for  the  arguments  leading  to  Eq.  (4)  to 
hold  /  c  must  be  larger  than  the  scales  of  surface  wrinkling  or  the  surface  contained  in  /  c3 
will  not  be  representative  of  all  scales  of  wrinkling.  At  the  same  time  /  c  cannot  be  too  large 
or  the  isotropic  character  will  be  lost.  Thus  /  c  =  /  =  e0  appears  to  be  the  obvious  choice. 
This  choice  for  /  c  implies  that  /  p  >  /  which  may  surprise  some  people.  Instantaneously 
the  flamelet  is  not  wrinkled  sufficiently  to  occupy  a  volume  comparable  to  the  entire 
turbulent  flame  brush.  Instead  it  occupies  at  any  instant  a  smaller  volume,  and  it  occupies 
the  entire  flame  brush  only  in  an  average  sense. 

An  alternative  line  of  argument  can  be  developed  which  leads  to  the  same  closure 
expression  as  above.  In  this  case  /  c//  F  is  set  equal  to  a  constant.  Thus  the  probability  of 
the  flamelet  being  in  /c3  for  a  given  <c>  but  different  /F  is  fixed.  lc  varies  with  /  F  and  the 


flamelet  does  not  fill  /c3.  Following  arguments  similar  to  those  of  [1]  one  concludes  that 
the  conditional  mean  flamelet  area  in  /  c3  is  given  as 
A;  =  (eo/Ei^Ao 

where  £j  =  t|/f,  £q  =  / ,  and  A0  is  defined  in  Fig.  1.  (Conditioning  is  on  the  presence  of  the 
flamelet  in  the  volume  defined  by  /  c3. )  Furthermore  A 0~/  c2  so  that 
A~(l  f/rj)0-2/  c2  and 
A//c3~(/  f/ri)0-2/,.-1 

which  leads  to  the  same  result  for  <cuf>  as  above  when  the  probability  Pc  is  incorporated. 
In  this  case  Pc  -  <c>(l-<c>)  l  cll  p- 

The  new  step  in  this  argument  relative  to  the  reasoning  of  [1]  is  the  scaling  A„  -  /c2. 
A0,  the  area  measured  for  e  2  e0,  is  the  mean  area  that  the  flamelet  would  have  while  in  /  c3 
if  the  flamelet  was  viewed  at  large  £  and  appeared  smooth,  i.  e.,  without  wrinkles.  In  this 
case  the  the  flamelet  would  be  planar  and  its  area  in  /  c3  scales  as  /  c2. 

In  summary  ,  the  final  expression  for  <%>  is 

<wf>  =  CRp0  <AYfUL>f  (/  (/  t\ )°-2l  f'1  <c>(l-<c>).  (2’) 

f  is  given  by 

f  =  [1  -  (1  -  At'l/4R/'3/4)exp(-At1,4R/'1/4U7<UL>f)]. 

Also  note  that  for  homogeneous,  isotropic  turbulence  /  /  rj  =  A^R/3'4  [1]. 

III.  Analyses  of  Normal  and  Oblique  Flames 

Normal  and  oblique  flames  are  analyzed  using  the  proposed  rate  model  as  presented 
in  (9).  Results  for  the  burning  velocity  are  compared  to  the  burning  velocity  model  of  [1]. 

For  a  normal  flame  assuming  high  Reynolds  number,  one-dimensional  flow  in  the 
mean  and  AYf  constant  (for  convenience),  one  can  develop  an  equation  for  <c>: 
p0utd<c>/dx  + 
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d/dx  (<pxc'u>  +  <cxp'u>  +  <Uxc'p>  +  <p’c'u>)  (6) 

=  <cof>/AYf. 

ut  is  the  mean  reactant  flow  velocity  approaching  the  flame  brush  or  time  mean  reaction 
zone  and  is  the  turbulent  burning  velocity.  Primes  denote  fluctuations  in  c  and  p,  while  U 
and  u  are  the  instantaneous  and  fluctuating  velocities  respectively.  If  this  equation  is 
integrated  from  -<*>  to  +°°  with  <c>  =  0  at  x  =  -  and  <c>  =  1  at  +oo,  one  obtains 

oo 

pout  =  J  <£Of>/AYf  dx, 

and  substituting  for  <cof> 

oo 

Pout  *  P0  {  cR<uL>f  (/cf/ti)D'2<c>(l-<c>)  /  p1  dx. 


It  is  implicit  in  the  formulation  presented  in  the  previous  section  that  most  of  the  spatial 
variation  in  <tof>  is  due  to  <c>(l-<c>),  i.e.,  it  is  in  the  probability  of  finding  a  flamelet  in 
/c3  at  x. 

Assume  that  <uL>f,  rj,  f  and  /  c  are  uniform  in  space  and  dependent  on  conditions  in 
the  reactants.  Also  note  that  the  fractal  modeling  in  [1]  gives  ut/<uL>f  =  (/  f/q)D'2.  It 
follows  that 


A  f. D-2  f  ..  .  . 

(— )  =  CR  ( — )  L  <c>(l-<c>)  dx 

T[  ■  q  J 


(7) 


For  this  relationship  to  be  valid  /  c  =  / ,  as  suggested  in  the  previous  section,  and 


oo 

CR  =  /  p/  J  <c>(l-<c>)  dx. 


(8) 


Clearly, 
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oo 

J  <c>(l-<c>)  dx 

-oo 

is  a  measure  of  the  turbulent  flame  brush  thickness.  Thus 

oo 

/  p/  J*  <c>(l -<c>)  dx 

— oo 

is  a  constant  for  a  given  geometry  and  initial  conditions.  The  application  of  the  present 
model  to  the  normal  flame  geometry  is  thus  seen  to  be  self-consistent  with  previous  work 
[1]. 


The  above  analysis  for  normal  flames  can  be  extended  to  oblique  flames  without 
great  difficulty  and  it  is  interesting  to  see  if  the  notion  of  a  turbulent  burning  velocity  may 
be  applied  to  oblique  flames. 

As  above,  assume  constant  density  and,  hence,  straight  streamlines.  We  are 
interested  in  cases  where  the  flame  brush  thickness  increases  in  the  streamwise  direction  as 
is  observed  experimentally,  and  we  assume  a  priori  that  the  mean  flame  brush  structure 
follows  similarity  at  least  locally.  We  assume  that  <c>  constant  surfaces  can  be 
approximated  locally  by  plane  surfaces  and,  to  first  order,  that  the  flame  brush  thickness 
grows  linearly  with,  distance  along  the  flame  brush.  Finally  we  assume  slow  growth  in 
flame  thickness  and  parabolic  flow. 

The  geometry  to  be  studied  is  shown  in  Fig  2.  Note  that  by  definition  the  £  axis 
lies  along  <c>  =  0.5,  and  T]  is  perpendicular  to  £,  while  the  x  and  y  axes  lie  along  and 
perpendicular  to  the  approach  flow,  respectively. 

A  control  volume  analysis  is  pursued  and  the  net  mean  flux  of  products  out  of  the 
control  volume  is  set  equal  to  the  mean  rate  of  production  in  the  volume.  For  a  rectangular 
volume  of  unit  depth  the  rate  of  product  formation  can  be  written  as 
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P  =  J  J  p<cof>dridC, 


r)R  and  %  define  the  edges  of  the  flame  brush,  and  p  is  the  ratio  of  the  amount  of  product 
mass  formed  per  unit  mass  of  fuel  reacted.  By  referring  to  the  above  model  for  <cof>  one 
notes  that,  if  <c>  satisfies  similarity,  that  is,  <c>  is  a  function  only  of  ri  /  (£  +  £g)  where 
is  a  suitably  chosen  constant,  and  e0/ei  is  constant,  the  integral  of  <cof>  over  rj  is  a  constant 
and  production  can  be  written  as 

P  =  8pl  , 

6) 


which  is  constant,  ie.,  independent  of  £. 


For  the  product  fluxes  into  and  out  of  the  control  volume  only  convection  is 
considered  for  faces  1  and  2,  and  turbulent  fluxes  across  the  other  two  faces  are  zero  since 
conditions  are  uniform  in  the  vicinity  of  both  faces.  Thus  only  convective  fluxes  need  be 
estimated  and  no  turbulence  modeling  is  needed  at  this  point.  The  neglect  of  turbulent 
fluxes  of  product  across  1  and  2  implies  a  boundary  layer  approxamation,  ie.,  parabolic 
flow  and  small  d 1  p/d£.  Product  flux  into  the  control  volume  is 
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Min  =  p<U>  cos0  [  /  F(l)  I£  +  tan(9  +  0p)  8],  where  similarity 
has  been  used  and 

^R 

Ic  =  J  <c>  drj  /  /  F 

The  angles  9  and  0P  are  defined  in  Fig.  2  and  /F(l)  denotes  the  flame  brush  thickness  at  1. 

For  the  flux  out  one  obtains 

Mout  =  p<U>  cos0  [/  p(2)  Ic  +  tan  0  5]. 

The  second  term  in  this  expression  is  the  flux  across  the  control  volume  face  which  is 
totally  in  products. 

Equating  the  difference  of  Moul  and  to  production  one  obtains 

5  pl^  =  cos0  p<U>  [  (/  p(2)  -  /  p(l))  Ic  +  (tan0  -  tan(0  +  0P))  5  ]. 

For  unit  Schmidt  number  and  <iOf>  an  even  function  of  q,  by  symmetry  0R  -  0  = 

0  +  0P  =  A0.  Then 

pIM  =  sin0  p<U>  [  1  +  (2IC  -  1)  tanA0  /  tan0  ],  (9) 

where  we  have  made  use  of  the  relationship  (/  p(2)  -  /  P(l))  =  25  tanA0.  Without  chemical 
reaction,  symmetry  requires  that  0  =  0  and  ^  =  1/2.  Here  we  look  for  a  solution  which 
gives  a  meaningful  definition  for  the  turbulent  burning  velocity.  One  can  obtain  such  a 
result  in  two  ways:  assuming  either  Ic  =  1/2  or  A0  =  0.  In  either  case  the  same  expression 
for  ut  is  obtained. 

put  =  sin0  pu  =  pIM  ,  (10) 

Thus  a  meaningful  burning  velocity,  defined  as  the  component  of  the  mean  approach  flow 
velocity  perpendicular  to  C,  is  obtained  and  is  the  integral  of  the  mean  product  production 

rate  across  the  flame  brush. 

The  control  volume  width,  5,  does  not  appear  in  the  final  expression  for  ut  as 
would  be  expected.  However,  implicit  in  the  analysis  is  the  assumption  that  similarity  is 
valid  for  a  control  volume  with  finite  5. 


Using  the  closure  model  presented  above  one  can  evaluate  Iw  and  finds  by 
substitution 

ut  =  CR  <uL>f  (f/  /n)D'2  pIr,  (11) 

where 

IR  =  J<c>(l-<c>)  drj  /  Ip- 

From  the  analysis  of  normal  flames  Cr  =  I/Ir  ,  and  one  sees  that  ut  as  defined  here  for 
oblique  flames  is  equal  to  that  for  normal  flames  and  that  utis  independent  of  £  provided 
similarity  holds  which  implies  nondecaying  turbulence. 

The  major  assumption  used  to  obtain  this  result  is  that  of  similarity.  This 
assumption  plays  a  major  role  in  our  evaluation  of  Min  and  M0lU.  It  should  be  noted,  on 
the  other  hand,  that  the  constancy  of  1^  is  inherent  in  the  modeling  of  <cof>,  provided  £0/ej 

is  constant  across  the  flame  brush,  and  does  not  require  an  additional  assumption.  The 
quantity  <c>(l-<c>)//  F  is  proportional  to  the  probability  of  finding  the  flamelet  at  a  point 
where  the  mean  degree  of  reaction  is  <c>,  and  thus  the  integral  IR  is  constant  because  the 
probability  of  finding  the  flamelet  somewhere  across  the  flame  brush  is  unity. 

The  condition  A0  =  0  corresponds  to  a  constant  flame  brush  thickness  which  is 
contrary  to  observations,  while  with  Ic  =  1/2  the  thickness  continues  to  increase  with 
increasing  If  <C0f>  is  an  even  function  of  rj,  by  symmetry,  we  expect  Ic,  which  is  the 

mean  value  of  <c>  when  averaged  across  the  flame  brush,  to  be  approximately  a  half. 
(N.B.:  our  numerical  calculations  give  Ic  =>  1/2.)  A  continually  growing  flame  brush  might 

be  thought  to  imply  an  increasing  combustion  rate  because  of  an  increase  in  flamelet  area. 
However,  as  noted,  our  rate  model  scales  <cof>  by  /  p,  and  I^,  is  independent  of  £.  This 

scaling  arises  via  Eq.  (4)  according  to  which  the  probability  Pc  of  finding  a  flamelet  in  a 
volume  /  c3  decreases  as  C,  increases.  A  physical  interpretation  of  this  behavior  is  that  the 
flame  brush  thickness  grows  with  increasing  C,  due  to  the  cumulative  effects  of  the  many, 
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random  upstream  displacements  of  the  flamelet  surface(s)  in  a  process  which  may  have 
some  likeness  to  a  random  walk.  These  displacements  cause  the  flame  brush  to  grow  but 
they  do  not  increase  flamelet  area  and  thereby  combustion  rate.  Flamelet  area  is  increased 
by  the  fractal  wrinkling  of  flamelet  surfaces  due  to  the  reactant  flow  turbulence.  This  area 
generation  is  dependent  on  /  f/rj  and  not  the  flame  brush  thickness. 
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IV.  Numerical  Analysis 


Results  of  calculations  using  the  proposed  rate  model  of  Eq.  (2')and  the  second 
order  transport  model  with  conditional  statistics  of  Chen,  Lumley  and  Gouldin  [4]  are  now 
presented  and  discussed.  The  transport  model  is  based  on  model  equations  for  conditioned 
moments,  where  conditioning  is  for  being  either  in  reactants,  c=0,  or  in  products,  c=l. 
The  instantaneous  progress  variable  acts  as  the  indicator  function. 

An  equation  for  the  mean  progress  variable,  for  negligible  density  change  and  high 
Reynolds  number,  can  be  written  [4]  as 


8<c>/3t  +  <U>«V<c>  =  V«{(<U>  -  <U>ri))<c>}  -  <o)f>. 

<U>ri)  is  the  conditioned  velocity  in  the  products,  and  the  first  term  on  the  right  hand  side 
of  the  equation  is  the  turbulent  flux  of  products.  This  term  would  be  given  by  -V»DtV<c> 
if  a  gradient  transport  model  were  used;  Dt  is  the  turbulent  diffusivity.  Instead  we  solve  a 
model  equation  for  <U>ri)  as  well  as  finding  <U>,  and  consequently,  gradient  transport 
modeling  is  avoided. 

These  model  equations  (  consult  [4]  for  more  detail  of  the  model )  are  combined 
with  the  new  rate  model,  and  calculations  are  performed  for  an  oblique  flame  in  decaying 
grid  turbulence.  The  flow  and  flame  configuration  are  shown  schematically  in  Fig.  3.  This 
is  the  same  configuration  as  considered  in  [4]  and  comparison  of  tfie  present  results  with 
those  in  [4]  can  be  made  by  the  reader  to  see  the  influence  of  changing  rate  models.  For 
calculations  the  flow  is  assumed  to  be  parabolic.  Upstream  boundary  conditions  and  model 
constants  for  the  calculations  are  given  in  Table  I;  note  that  D  =  2.4. 

As  a  test  of  the  turbulence  model  for  nonreacting  flow,  the  decay  of  grid  turbulence 
was  calculated,  and  the  results  for  <u2>,  /  and  R/  compared  to  the  experimental  data 

ofWarhaft  [6].  Comparisons  are  made  in  Fig.  4  were  it  can  be  seen  that  model  predictions 
for  <u2>  and  R<  are  very  good. 
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If  the  assumptions  of  the  analysis  of  the  previous  section  are  valid  one  would 
expect  that  the  turbulent  burning  velocity  defined  as  the  component  of  the  mean  reactant 
flow  velocity  normal  to  the  <c>  =  0.5  surface  would  equal  that  predicted  by  the  fractal 
model  of  [1];  see  Eq.  (11).  In  Fig.  5  calculation  results  for  ut/u0  are  compared  to  those 
predicted  by  the  burning  velocity  model.  For  this  comparison  CR  is  chosen  to  give  the  best 
overall  agreement;  CR  =  4.  The  comparisons  shown  in  Fig.  5  for  a  broad  range  of 
conditions  are  quite  good  both  with  regard  to  magnitude  and  variation  of  ut/u0  with  x/M. 

We  attribute  the  small  observed  variations,  which  are  present,  between  burning  velocity 
model  and  calculation  results  to  a  departure  from  flow  similarity  which  is  assumed  in  the 
analysis  of  the  last  section. 

A  lack  of  similarity  can  be  seen  in  the  variation  of  lateral  <c>  profiles  with  x/M,  and 
several  typical  profiles  are  shown  in  Fig.  6.  By  similarity  one  would  expect  to  see  /  p  grow 
linearly  with  x/M.  Such  behavior  is  not  found  in  the  numerical  results.  Instead  /  F  is  found 

to  grow  very  rapidly  at  small  x/M  and  then  to  grow  much  more  slowly  at  larger  x/M;  see 
Fig.  7,  which  presents  /  F/  x  versus  x  /  M. 

In  performing  the  calculations  for  this  work  some  results  were  also  obtained  with  a 
gradient  transport  model  for  the  <c>  equation  while  a  full  second  order  model  was  retained 
for  the  turbulent  Reynold  stresses.  For  these  calculations  Dt  =  Cdk2/e,  where  k  is  the 
turbulence  kinetic  energy  and  e  is  the  turbulence  energy  dissipation  rate;  the  constant  Cd  is 

set  equal  to  0.09  [7].  We  find  that  results  for  burning  velocity  are  similar  for  the  two 
transport  models,  Fig.  8,  but  that  there  is  a  marked  difference  in  /  F,  Fig.  9.  From  Fig. 

10,  it  can  be  seen  that  gradient  transport  greatly  under  predicts  turbulent  fluxes  in 
comparison  with  the  second  order  model. 

Counter-gradient  diffusion  is  not  observed.  Libby,  et  al.  [8]  have  shown  that  the 
interaction  between  the  mean  pressure  gradient  and  density  fluctuations  plays  a  crucial  role 
in  causing  counter-gradient  diffusion.  Pressure  gradient  fluctuations  may  also  contribute  to 


this  phenomenon.  Thus  it  is  not  surprising  that  when  combustion  with  negligible  density 
change  is  considered  counter-gradient  diffusion  is  not  observed. 

That  both  transport  models  give  essentially  the  same  burning  velocity  may  at  first  be 
surprising.  The  reason  for  this  behavior  is  as  follows.  From  Eqs.  9  and  10  we  see  that 
with  similarity  put  is  given  by  pl^  for  Ic  =  1/2.  As  discussed  above,  with  our  rate  model  Iu 

is  dependent  only  on  <p0AYf  uL>  (f l  /rj)0-2  and  hence  independent  of  the  species  transport 
model.  This  result  follows  from  our  expression  for  Pc  and  its  scaling  by  /  p'1.  Then  since 
both  transport  models  give  Ic  ~  1/2,  they  give  similar  ut  values. 

Conditional  and  unconditional  moments  of  velocity  are  presented  in  Figs.  1 1  and 
12.  It  can  be  seen  that  unconditional  mean  velocities  in  the  transverse  (y)  direction  remain 
close  to  zero.  They  should  equal  zero  in  the  absence  of  density  change.  Conditional  mean 
velocities  (  Fig.  11  )  for  the  burned  and  unbumed  zones  are  of  equal  magnitude  but 
opposite  sign,  and  their  difference  decreases  in  the  flow  direction  as  the  turbulence  decays. 
The  relative  magnitude  and  sign  of  the  conditional  transverse  velocities  must  undergo 
significant  change  in  the  presence  of  combustion  with  heat  release  and  density  change  if 
counter-gradient  diffusion  is  to  be  observed.  Transverse  turbulent  normal  stress  profiles 
are  presented  in  Fig.  12.  The  unconditioned  stresses  are  nearly  uniform  as  is  required  in  the 
absence  of  density  change,  while  conditional  stresses  vary  significantly  across  the  flame 
brush. 

From  the  above  presentation  and  the  figures  it  is  clear  that  the  new  rate  model  gives 
consistent  and  plausible  results.  Predicted  burning  velocities  are  in  the  expected  range. 
However  exact  correspondence  to  the  burning  velocity  model  of  Section  III  is  not  achieved. 
This  lack  of  agreement  is  most  likely  the  result  of  the  departure  from  flow  similarity  in  the 
numerical  results  which  in  turn  we  attribute  to  the  decay  of  turbulence. 

In  the  model  calculations  a  constant  density  has  been  assumed  and  no  counter¬ 
gradient  diffusion  is  observed.  At  present  the  conditional  modeling  approach  has  not  been 


adapted  to  account  for  density  changes  across  the  flamelet.  Work  to  allow  for  density 
change  is  needed. 
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V.  Summary  and  Closing  Comments 

A  closure  model  for  mean  fuel  consumption  rate  in  premixed  turbulent  flames  has 
been  proposed  and  applied  to  several  test  problems.  The  model  expresses  the  mean  rate  in 
terms  of  the  product  of  the  mean  fuel  consumption  rate  per  unit  flamelet  surface  area  and  of 
an  expression  for  the  mean  flamelet  area.  The  latter  expression  is  obtained  under  the 
assumption  that  flamelet  surface  geometry  can  be  represented  by  fractal  surfaces.  In 
addition  the  model  uses  an  empirical  expression  from  [5]  for  the  probability  of  finding  the 
flamelet  in  a  volume  /c3. 

The  new  parameter  appearing  with  the  model  is  the  fractal  dimension.  Based  on  the 
results  reported  in  [1]  and  related  experimental  [9,10]  and  theoretical  results  [1 1,12,13]  we 
recommend  that  D=2.37,  although  for  the  present  calculations  we  have  used  a  value  of  2.4. 

Test  problem  results  are  mutually  consistant  and  pi;  asible.  The  rate  model  is  easy 
to  use  and  gives  burning  velocity  values  which  are  equal  to  the  predictions  of  a  new 
burning  velocity  model  [1]  which  also  uses  fractal  concepts. 

The  oblique  flame  analysis  suggests  that  a  meaningful  burning  velocity  can  be 
defined  for  oblique  flames,  and  that  ut  values  can  be  related  to  the  traditional  definition  of 

the  burning  velocity  of  normal  flames.  However,  the  oblique  flame  analysis  rests  on  the 
assumption  of  similarity,  and  the  expression  for  burning  velocity  requires  a  supplemental 
condition  --  either  A0  =  0  or  Ic  =  1/2.  The  numerical  calculations  of  Section  IV  show  that 

departures  from  similarity  can  affect  the  burning  velocity  variation.  In  the  absence  of 
similarity  it  may  be  possible  to  develop  an  alternative,  useful  definition  for  the  burning 
velocity.  Until  such  a  definition  has  been  established,  turbulent  burning  velocity  data 
should  be  interpreted  with  caution. 

The  empirical  basis  for  the  relationship  Pc~<c>(l-<c>)  can  and  should  be 

made  more  substantial.  For  example  the  passage  time  model  of  [14]  suggests  P~<o(l- 
<c>)  but  not  the  normalization  by  /F'*.  Furthermore  in  [5]  the  turbulence  Reynolds  number 


is  low  and  multiple  flamelet  crossings  of  lines  parallel  to  the  r|  axis  are  rare  events.  It  is 
reasonable  to  expect  that  future  work  will  lead  to  an  improved  expression  of  Pc. 
Developments  of  the  passage  time  analyses  of  Bray,  Libby  and  Moss  [14,  15]  may  help  in 
this  respect.  The  possible  occurrence  of  large  scale  doubling  as  described  in  Section  II 
needs  further  study.  If  such  events  are  important  it  may  be  necessary  to  change  our  scaling 
of  /c  and  our  expression  for  Pc.  To  avoid  confusion  it  should  be  noted  that  the  neglect  of 

large  scale  doubling  back  does  not  preclude  doubling  back  on  scales  up  to  the  order  of  the 
outer  cutoff.  Thus  along  lines  parallel  to  the  r\  axis  multiple  flamelet  crossings  at  any 
instant  will  occur  over  a  length  of  order  of  the  outer  cutoff.  The  region  of  these  crossings 
will  fluctuate  in  time  and  thus  form  a  flame  brush  which  is  thicker  than  the  outer  cutoff. 

The  ratio  r|/f  has  been  used  to  represent  the  inner  cutoff.  There  is  some  question  as 
to  the  appropriateness  of  this  representation  [1,16].  Our  expression  for  rj/f  has  been  taken 
from  [1]  and  the  present  work  provides  no  test  of  its  validity.  Should  a  better 
representation  of  e*  be  developed  it  is  easily  incorporated  in  the  model. 

The  fractal  character  of  flamelets  is  open  to  direct  experimental  verification  by 
methods  such  as  those  described  by  Sreenivasan  and  Maneveau  [9].  Information  on  the 
fractal  dimension  and  the  cutoffs  can  also  be  obtained.  Results  available  to  date  are 
preliminary  [16,17],  and  they  are  for  low  turbulence  Reynolds  numbers.  While  they  show 
fractal  behavior,  inferred  D  values  are  low  --  between  2.1  and  2.2.  These  low  values  may 
be  related  to  the  low  turbulence  Reynolds  number  of  the  flow  and  to  the  lack  of  a  developed 
inertial  subrange  in  the  velocity  fluctuation  power  spectrum.  Thus  measurements  at  high 
turbulence  Reynolds  numbers  are  required. 
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Variation  with  measurement  scale  of  measured  area  in  a  cubic  volume  of 
side  L  for  a  surface  exhibiting  fractal  character,  a)  No  inner  or  outer 
cutoffs,  b)  Inner  and  outer  cutoffs  to  fractal  behavior. 

Schematic  diagram  of  flame  brush  showing  coordinate  systems. 

Sketch  of  flow  configuration  assumed  for  calculation. 

Comparisons  of  model  predictions  for  (a)  <uu>/<U2>  and  (b)  R/ with 
Warhaft  data  [6]. 

Comparisons  of  burning  velocities  obtained  with  new  rate  model  and 
predictions  from  model  of  Reference  1:  a)  varying  R/  (by  varying  v),  b) 
varying  unstrained  laminar  burning  velocities,  u0  =  2uT>  ur  u/2  and  1^  =  1 

cm/s,  and  c)  varying  initial  turbulence  integral  scale,  lr  =  7.27mm  (note  v  is 
adjusted  to  give  constant  Rt). 

<0  profiles  at  three  locations  for  standard  conditions.  X/M  =  50,100  and 

200. 

Flame  brush  thickness  for  different  R/_  see  Fig.  5a  for  burning  velocity 
results. 

Burning  velocity  predictions  from  second-order  model  and  from  gradient 
transport  model,  standard  conditions.  Burning  velocity  model  predictions 
are  also  shown. 

Flame  brush  thickness  predictions  for  second-order  model  and  gradient 
transport  model,  standard  conditions. 

Turbulent  fluxes  for  standard  conditions  a)  firm  conditional  statistics  and 
b)  from  gradient  transport 

Second  order  model  predictions  of  transverse  mean  velocity,  both 
conditioned  and  unconditioned.  Note  that  <vc'>  =  <c>(l-<c>)  (<V>0)- 
<V>(°)).  x/M  =  50.0. 

Second  order  model  predictions  of  transverse  turbulent  normal  stresses, 
both  conditioned  and  unconditioned.  x/M  =  50.0. 
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TABLE  I 

Standard  Inlet  Conditions:  <u>  =  7  m/s 

ul<  U>  =  0.07 
u  =  49  cm/s 
lT  =  7.27  mm* 

Uj  =1  cm/s  (laminar  flame  speed)* 

v  =  1.85X10-5  m2/s 
R/=  229 

Model  Constants:  CR  =  4.0 

D  *  2.4 
At=  0.37  [1] 

Cd  =  0.09  [7] 

For  other  constants  see  [4] 

*  r  denotes  standard  reference  value. 
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Fig.  1.  Variation  with  measurement  scale  of  measured  area  in  a  cubic  volume  of  side  L 
for  a  surface  exhibiting  fractal  character,  a)  No  inner  cm- outer  cutoff,  b)  Inner 
and  outer  cutoffs  to  fractal  behavior. 
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An  Interpretation  of  Jet  Mixing  Using  Fractals 
F.  C.  Gouldin 

Mechanical  and  Aerospace  Engineering 
Cornell  University 
Ithaca,  New  York  14853 

ABSTRACT 


An  analysis  of  axisymmetric  jet  mixing  is  reported  in  which  the  jet  structure  is 
modeled  as  a  set  of  constant  concentration  surfaces.  Stationary  turbulence  is  assumed,  and 
a  balance  equation  for  the  ensemble  mean  flux  of  jet  fluid  across  constant  jet  fluid 
concentration  surfaces  is  developed.  Terms  in  the  equation  are  modeled  1)  by  assuming 
that  the  geometry  of  constant  concentration  surfaces  can  be  representerd  by  fractal  surfaces 
and  2)  by  estimating  the  mean  jet  fluid  mass  flux  per  unit  area  of  a  constant  concentration 
surface  as  the  Kolmogorov  velocity.  Flow  similarity  is  assumed.  It  is  found  that  for 
Reynolds  number  independence  the  fractal  dimension  is  given  by 


D  =  2  + 


pi  +  2 


In  addition  an  expression  for  the  centerline  variation  of  mean  jet  fliud  mass  fraction  is 
obtained;  the  expression  is  of  similarity  form  as  expected. 
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Introduction 


In  [  1]  Mandelbrot  proposes  that  constant  scalar  property  surfaces  in  homogeneous, 
isotropic  turbulence  are  fractal  surfaces  and  suggests  two  values  for  the  fractal  dimension 
depending  on  the  spatial  variance  of  the  scalar  field.  D  =  2  1/2  for  Gauss-Burger 
turbulence  and  2  2/3  for  Gauss-Kolmogorov  turbulence.  Lovejoy  for  clouds  [2]  and 
Sreenivasan  and  Maneveau  [3]  for  shear  flows  obtain  data  suggesting  fractal  dimensions 
between  approximately  2.32  and  2.4  for  constant  property  surfaces.  These  values  are 
noticeably  lower  than  the  two  suggested  by  Mandelbrot  Hentschel  and  Procaccia  [4]  and 
more  recently  Kingdon  and  Ball  [5]  present  analyses  of  cloud  dispersion  in  homogeneous, 
isotopic  turbulence.  In  these  analyses  molecular  transport  is  ignored,  high  turbulence 
Reynolds  number  is  assumed  and  the  dispersion  of  one  species  into  another  is  modeled  by 
considering  the  quantity  Q(J,  1 1  /  ',0),  which  is  the  probability  that  a  pair  of  diffusing 
particles  of  fluid  separated  by  / '  at  t  =  0  will  be  a  vector  distance  l  apart  at  t  Q  in  turn  is 
related  to  the  structure  function  S(/  ,t)  at  a  point  r  for  a  scalar  0: 

S(/,t)  =  —  J  dr  <|9(r  +  / ,  t)  -  0(r,  t)|2> 

For  small  \l  |  =  /,  S (/ ,  t)  scales  as  a  power  of  / ,  and  the  exponent  gives  the  fractal 
dimension. 

Both  Hentschel  and  Procaccia  and  Kingdon  and  Ball  develop  models  for  Q  and 
suggest  a  form  for  S (/  ,t)  as  /  ->  0. 

lint/  _*o  S(/ )  -  / 8 

with  B  constant.  Both  arguments  contain  heuristic  elements.  There  is  a  fundamental 
difference  between  the  two  sets  of  authors  in  their  modeling  of  Q  and  in  their  expressions 
for  lim/  S(/ ).  Yet  their  final  expressions  for  D  are  the  same: 
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D  =  2  + 


\x  +  2 

~T~ 


(i  is  the  intermittency  exponent. 

In  this  paper  a  semi-empirical,  physically  oriented  model  of  mixing  in  jet  flows  is 
developed  which  assumes  the  instantaneous  jet  structure  can  be  represented  by  an  ensemble 
of  fractal,  constant  property  surfaces.  The  model  gives  an  expression  for  mean  jet  fluid 
concentration  on  the  jet  centerline  which  contains  the  jet  Reynolds  number  raised  to  a 
power.  For  there  to  be  no  Reynolds  number  dependency  as  required  for  similarity,  it  is 
necessary  for  the  exponent  to  equal  zero.  As  a  result,  one  finds  D  =  2  +  ((2+n)/6)  as 
obtained  both  by  Hentschel  and  Procaccia  and  by  Kingdon  and  Ball.  This  expression  for  D 
is  the  major  finding  to  be  reported.  While  it  is  identical  to  that  of  others  it  is  obtained  by 
very  different  reasoning  and  for  a  different  flow,  a  free  shear  flow. 

Jet  Mixing  Model 

For  analysis  the  following  major  assumptions  are  made.  The  flow  is  an 
axisymmetric  jet  into  a  stagnant  fluid  with  uniform  and  constant  density.  The  turbulence  is 
assumed  to  be  stationary  and  jet  Reynolds  number  is  large  such  that  the  flow  is  Reynolds 
number  independent.  The  distributions  of  mean  composition  and  mean  velocity  follow 
similarity.  It  is  further  assumed  that  surfaces  of  constant  jet  fluid  mass  fraction  may  be 
represented  by  fractal  surfaces  for  the  purpose  of  estimating  the  ensemble  mean  area  of 
such  surfaces. 

With  similarity  the  mean  jet  structure  is  given  by  the  following  expressions  [6]. 


<U>  =  u0  f(r/r1/2). 


(1) 
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Uj  /u0  =  Cu  x/d. 

(2) 

<Z>  =  Zq  g(r/r1/2). 

(3) 

1/  Zo  =  Cz  x/d. 

(4) 

The  quantities  in  the  above  equations  are  defined  as  follows:  x  is  axial  distance;  r  is  radial 
distance.  U  is  the  axial  velocity,  Uj  its  momentum  flux  weighted  average  at  the  jet  origin 
(x=0),  and  u0  its  mean  value  on  the  jet  axis.  Z  is  the  jet  fluid  mass  fraction  and  Zq  its  mean 
value  on  the  axis,  d  is  the  initial  jet  diameter,  and  r1/2  is  the  jet  half  radius  Cu  and  C,  are 
empirical  constants.  Finally  <  >  denotes  ensemble  mean.  (4)  gives  the  centerline  decay  of 
the  mean  jet  fluid  mass  fraction.  The  modeling  presented  here  will  lead  to  an  expression 
similar  to  (4)  with  a  term  R/*0*,  where  R;  is  a  turbulence  Reynolds  number  based  on  the 

integral  scale,  /,  and  the  root  mean  square  of  velocity  fluctuations.  For  similarity  h(D)  *  0 
and  an  expression  for  D  is  obtained. 


From  experiment  it  is  found  that  f  and  g  can  be  expressed  to  reasonable  accuracy  by 


exp(-(ar/rl/2)2),  (5) 

where  a  is  an  empirical  constant. 

The  instantaneous  jet  structure  can  be  viewed  as  a  set  of  constant  Z  surfaces  and  it 
is  hypothesized  that  these  surfaces  exhibit  fractal  character  over  a  range  of  length  scales 
associated  with  the  scales  of  the  turbulent  velocity  fluctuations,  namely  rj  ( the  Kolmogorov 
scale  )  to  / .  In  this  work  a  fractal  description  is  assumed  to  obtain  an  expression  for  mean 
surface  area.  The  success  of  the  resultant  model  supports  the  fractal  hypothesis  and  a 
certain  expression  for  D  but  does  not  prove  the  fractal  hypothesis. 


B5 


For  a  surface,  S,  enclosing  a  volume,  V,  an  integral  equation  for  jet  fluid 
conservation  may  be  written  as 

-^JpZdT  =  j”‘”di  +  JPZqb*ds  =  -J[pZ(q*qb)+jz]»  ds. 

V  v  ^  s  s 

This  is  a  general  equation  for  a  surface  moving  with  an  arbitrary  local  velocity,  qb.  q  is  the 
fluid  velocity,  and  jz  is  the  diffusive  flux  of  jet  fluid.  Let  the  surface  be  defined  by  a 
constant  value  for  Z,  and  take  the  ensemble  average  of  this  equation.  For  stationary 
turbulence  the  mean  of  the  left  hand  side  must  be  zero,  and  one  obtains 

<JpZ  (q  -  qb)  •  ds  >  +<Jjz»ds>  =0  (6) 

s  s 


For  jet  flow  the  surface  integrals  may  be  written  as  the  sum  of  two  types  of  terms  -  a  flux 
of  fluid  from  the  jet  and  a  surface  integral. 

\T  =  Jt(d/2)2  pV.  =  <|[pZ  (q  -  qb)  +  jz]  •  ds  >.  (7) 

s 

Vj  is  the  initial  volume  flux  weighted  average  jet  velocity.  The  areas  S  in  (6)  and  (7)  differ 
by  the  initial  cross-sectional  area  of  the  jet.  It  is  assumed  that  Z  constant  surfaces  are 
attached  to  the  lip  of  the  jet  and  that  these  surfaces  are  simply  connected.  The  first 
assumption  is  valid,  and  if  the  second  assumption  is  relaxed  it  seems  reasonable  to  argue 
that  the  average  result,  (7),  is  not  affected.  The  right  hand  side  of  (7) 
is  now  evaluated  approximately. 


Consider  the  segment  of  Z  constant  surface  between  x  and  x+dx.  We  approximate 
the  contribution  from  the  surface  in  this  region  to  the  right  hand  side  of  (7)  by  a  mean  flux 
times  a  mean  surface  area.  (7)  is  then  estimated  by  integrating  these  contributions  over  x. 
The  mean  surface  area  is  estimated  with  the  help  of  the  fractal  hypothesis. 

As  noted  above  we  assume  fractal  behavior  of  Z  surfaces  for  a  range  of 
length  scales  between  an  inner  cutoff  (£j)  and  an  outer  cutoff  (e0).  Measurements  of  area  at 

scales  outside  the  range  of  fractal  behavior  are  assumed  to  give  one  of  two  constant  values 
corresponding  to  large  or  small  measurement  scale.  In  [7],  it  is  argued  that  Aj,  the  area  of  a 

fractal  surface  measured  at  scales  below  the  inner  cutoff,  is  the  actual  ensemble  mean 
surface  area  of  a  constant  property  surface,  i.e.,  Az  =  A,.  Also  from  [7]  one  may  write  for 

the  ratio  of  the  two  surface  areas  measured  beyond  the  cutoffs 

Ai/A0  =  (e0/ei)D-2  =  (//T1)D-2t  (8) 

where  l  and  r\  are  used  as  the  outer  and  inner  cutoffs.  The  choice  of  rj  as  the  inner  cutoff 
implies  a  unit  Schmidt  number.  For  homogeneous,  isotopic  turbulence  /  /rj  =  At1/4  Rj 3/4 
[7]  with  At  a  constant,  and  substitution  gives 

Aj/A0  =  (At1/4  R/ 3/4  )D'2  .  (9) 

A0  is  the  surface  area  measured  for  e  >  / .  By  assumption  A0  for  a  Z  surface  is 
equated  to  the  area  of  the  <Z>  constant  surface  for  which  <Z>  equals  Z.  Sketches  of  the 
intersections  of  a  <Z>  constant  surface  and  of  realizations  of  a  Z  constant  surface  with  a 
plane  containing  the  jet  axis  are  depicted  in  Fig.  1.  Note  in  Fig.  1,  L  is  defined  as  the  axial 
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location  where  the  <Z>  constant  surface  crosses  the  x  axis.  Eq.  9  is  now  applied  to  obtain 
an  expression  for  dAz,  the  mean  Z  surface  area  between  x  and  x+dx. 

The  area  of  the  <Z>  surface  in  the  region  x  to  x+dx  is 

dA<Z>  =  2jrr<Z>dx« 

and  therefore  if  dAz /  dA<Z>  =  A^/Aq,  one  obtains 

d AZ  =  2nr<Z>  (Ati'4  R,  3/4)D‘2dx  (10) 

as  the  ensemble  average  of  the  Z  surface  area  in  the  region  x,  x  +  dx.  r<z>  is  the  radial 
location  of  the  <Z>  surface  at  x. 

The  arguments  leading  through  (8)  to  (10)  may  be  made  graphically.  For  a  surface 
exhibiting  fractal  behavior  between  length  scales  E;  to  e0,  the  measured  surface  area  is 

assumed  to  vary  with  measurement  scale  as  shown  in  Fig.  2.  Fractal  behavior  is  given  by 
the  power  law  dependency  of  A  upon  e.  For  e— >0,  A  must  approach  a  constant  value,  the 
true  value,  if  the  surface  has  a  minimum  scale  of  wrinkling  as  assumed.  For  a  Z  surface 
molecular  diffusion  is  expected  to  limit  the  minimum  scale  of  wrinkling.  This  limiting 
behavior  at  small  e,  is  represented  by  the  horizontal  line  in  Fig.  2  at  A  =  At.  We  define  £j 
by  Aj/Lj  -  (£j/L)2  *  D/L  and  assume  on  physical  grounds  that  £4  =  r|.  ( N.B.:  since  the  ratio 
A/A0  is  required,  the  constant  of  proportionality  is  unimportant  in  this  application. ) 

The  behavior  of  A  versus  e  depicted  in  Fig.  2  for  large  £  implies  that  the  surface 
appears  smooth  and  planar  at  large  length  scale.  The  area  of  this  planar  surface  is  A0,  and 
the  outer  cutoff  is  defined  by  AjV  -  (£0/L)2  D/L.  We  assume  that  £0  =  / .  For  this  planar 
assumption  to  be  valid  it  is  expected  that  £0  must  be  much  less  than  the  radius  of  curvature 
of  the  <  Z  >  surface.  This  limiting  behavior  for  large  £  is  reasonable  for  shear  flows  but 
may  not  be  reasonable  for  homogeneous,  isotopic  turbulence. 


Physically  it  is  expected  that  the  Z  surface  will  spend  most  of  its  time  in  the  vicinity 
of  the  <  Z  >  surface.  Hence  A0  is  identified  with  the  latter  surface. 

Finally,  it  is  assumed  that  (10)  refers  to  the  average  surface  area  between  x  and 
x+dx.  The  actual  surface  area  should  very  from  realization  to  realization,  and  the  surface 
for  each  such  realization  is  expected  to  exhibit  fractal  character.  Thus  the  identification  of 
(10)  with  mean  values  is  nontrivial  and  open  to  question.  A  similar  assumption  is  made  in 
[7]  with  success.  A  fuller  justification  of  the  assumption  is  left  to  future  work. 

The  flux  of  jet  fluid  across  a  Z  surface  is  the  result  of  convection  (q  -  qb)  and 
diffusion  (jz).  To  estimate  the  mean  of  these  two  fluxes  we  write  ud  =  D/"nd,  where  r|d  is 

an  effective  length  scale  defined  by  the  above  expression,  and  the  total  mean  mass  flux 
across  the  surface  per  unit  area  is  pZud.  (N.B.:  D  here  is  the  diffusion  coefficient  There 

should  be  no  confusion  with  the  fractal  dimension  since  the  latter  appears  as  an  exponent) 
Gibson  [8]  has  considered  the  motion  of  a  constant  property  surface  in  turbulent  flow. 
From  his  analysis 

(q  -  qb)  =  -D(V2Z/ 1 VZ  |  )(VZ / 1 VZ I ) 


in  the  present  context,  while  the  diffusive  flux  is  given  by 


jz  s  -  pZDVlnZ, 


for  Fick's  law  diffusion.  Thus  the  appropriate  length  scale  for  ud  is  related  to 

<  |  (V2Z/ 1  VZ !  )(VZ / 1  VZ  I )  +  VlnZ  I  >  ‘. 


One  expects  intuitively  that  the  largest  contributions  to  rjd  come  from  regions  in  the 
flow  where  Z  is  varying  rapidly  with  distance  and  hence  both  V2Z  and  VlnZ  are  large. 


These  regions  contain  the  smallest  turbulent  eddies  and  thus  an  appropriate  estimate  for  r|d 
would  seem  to  be  a  quantity  proportional  to  q  for  a  unit  Schmidt  number. 

Assuming  unit  Schmidt  number,  we  replace  D  by  v  and  estimate  rid  by  (1/Cd)ri, 
which  gives  ud  =  Cdv,  where  v  is  the  Kolmogorov  velocity  scale  and  Cd  is  a  model 
constant.  Furthermore,  for  homogeneous,  isotopic  turbulence  v  =  At1/4R/'1/4  u\  We 
will  use  this  expression  for  v  arguing  that  since  the  Z  surface  is  embedded  in  turbulent  fluid 
and  thus  u'  is  conditioned  on  being  in  turbulent  fluid  the  use  of  an  expression  for 
homogeneous,  isotopic  turbulence  is  reasonable.  Further  we  assume  that  u’  can  be 
expressed  in  terms  of  u0;  u'  =  (u7uo)u0  with  u'/u0  being  constant  over  the  jet. 
Conditioned  velocity  measurements  in  jets  show  that  u7u0,  while  not  constant  in  a  jet,  does 

not  vary  gready,  and  thus  the  above  assumption  is  reasonable  as  a  first  approximation.  An 
expression  for  ud  can  now  be  written. 

ud  =  CdAt1/4  Rj"1/4(u7u0)u0  .  (11) 

According  to  the  Gibson  expression,  q  -  qb  is  undefined  or  infinite  at  critical 
points,  i.e.,  at  points  where  VZ  =  0.  On  physical  grounds  one  expects  q  -  qb  at  critical 
points  to  be  infinite  and  to  be  large  in  the  region  of  a  critical  point.  Thus  it  seems 
reasonable  to  argue  that  surfaces  will  move  away  or  through  and  then  away  from  critical 
points  very  rapidly.  Hence  the  contributions  of  such  points  to  <q  -  qb>  are  expected  to  be 

negligible.  Alternatively,  one  can  argue  that  the  ensemble  average  of  the  integral  of  (q- 
qb)«ds  over  the  area  of  a  constant  property  surface  must  be  finite  even  if  q  -  qb  at  singular 
points  is  not.  Then  <q  -  qb>  is  defined  by  the  ratio  of  the  ensemble  average  integral  and 
the  ensemble  average  area.  The  second  agreement,  which  complements  the  first,  is  based 
on  the  physically  reasonable  assumption  that  the  average  convective  flux  of  fluid  across  a 


constant  property  surface  is  finite. 


By  combining  (10)  and  (11),  and  integrating  along  the  x-axis  one  obtains  an 
expression  for  Mj. 

L 

Mj  =  pZ  !  CdAt»'«  V1/4  (U7u0)u0  2rtr<Z>(Al1/4  R/3/4)D-2dx  .  (12) 

0 

To  proceed  further,  relationships  for  r<Z>  and  u0  as  functions  of  x  are  required. 

These  functions  depend  on  the  spreading  rate  of  the  jet,  and  the  basis  for  a  theoretical 
development  to  obtain  them  is  not  contained  in  the  fractal  picture  being  explored.  As  an 
alternative  to  such  a  development  empirical  expressions  based  on  similarity  arc  used  for  the 
needed  relationships.  Having  taken  such  a  step  one  cannot  claim  that  the  resulting 
expression  for  centerline  concentration  is  a  predictive  one.  The  exercise  does,  however, 
allow  one  to  determine  whether  or  not  the  fractal  picture  as  developed  gives  results 
consistent  with  experiment,  thereby  providing  a  test  of  the  applicability  of  fractal  concepts 
to  jet  mixing.  In  addition,  a  value  for  D  is  obtained. 

From  similarity  and  experiment  r1/2/d  =  (1/2)LZ  x/d  with  Lz  a  constant  [6].  With 
this  expression  and  (3)  and  (5)  one  can  write  for  r^,: 

r<z>  -  (d/2)  Lz  a'l/In (Ux)  x/d,  (13) 

where  use  has  been  made  of  the  fact  that  Zq  falls  to  <Z>  at  x  =  L.  (2),  (12)  and  (13)  can  be 
combined  to  give 

Vj  =  (Lj/Cua)  CdAt<D-i)'4R,(3/4(i>2)  -1/4)  (u7u0)Uj  • 

l  _ 

Z(L/d)jVln(l/Si)  d$,  (14) 

o 

where  %  =  x/L.  To  obtain  (14)  use  has  been  made  of  the  fact  that  Rt ,  by  similarity,  is 
independent  of  x.  The  integral  in  (14)  is  (rc/2)|/2. 


Uj  and  Vj  are  momentum  and  volume  average  velocities  and  their  ratio  will  very 
with  initial  condition;  let  Vj/Uj  =  b.  L  is  the  axial  distance  at  which  the  mean  jet  fluid  mass 
fraction  on  the  centerline  equals  Z.  Therefore  (14)  and  (4)  are  identical  provided 

h(D)  =  (3/4(D-2)-l/4)  =  0 


and 


(7C/2)‘/2(LJ/C»ab)CdAl(D-1V4(u7uo)  =  Cz. 

The  second  condition  is  reasonable.  (N.B.:  Cz  varies  with  initial  conditions  as  does  b 
[6].)  Substitution  with  reasonable  values  [6,  9]  into  the  above  gives  Cd  *  0.38.  The  first 
condition  is  an  equation  for  D. 

The  analysis  to  this  point  has  not  considered  the  intermittency  of  viscous 
dissipation.  The  effect  of  this  intermittency  is  to  alter  the  scaling  of  ud  and  of  rj//. 

Adjustments  in  the  analysis  are  straight  forward.  With  intermittency  the  dissipative  length 
scale  ld  replaces  as  the  inner  cutoff  and  from  Frisch,  et  al.  [10]  /d~  / 
(N.B.:Frisch,  et  al.  introduce  their  own  fractal  dimension,  De,  such  that  De  =  3  -  [1.  Here 

the  more  common  jj.  is  used  in  part  to  avoid  confusion  with  D  which  is  related  to  but  not 
equal  to  De.)  ud  is  scaled  as  Ed1/2  where  Ed  is  "the  kinetic  energy  per  unit  mass  on  scales" 
-  /  d  [10].  Again  from  [10]  it  can  be  shown  that  this  scaling  gives  ud  =  u’R/1/2  '3/(4*^. 
Thus  with  intermittency  the  Reynolds  dependency  in  (14)  becomes 

R/(3/(4-u)HD-2]R;l/2-3/(4-^)  ? 


and  for  Reynolds  number  independence 


D  =  2  +  (|i  +  2)16. 


(15) 
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This  result  is  identical  to  that  obtained  by  Hentschel  and  Procaccia  [4]  and  by 
Kingdon  and  Ball  [5] .  Our  result  for  jj.  =  0  is  not  equivalent  to  Mandelbrot's  result  (  D  =  2 
2/3)  for  Kolmogorov  turbulence  [1]  even  though  in  both  cases  the  variance  scales  as 
separation  distance  to  the  2/3  power.  In  his  development,  Mandelbrot  assumes  a  Guassian 
random  field  ;  no  such  assumption  is  made  in  the  present  development.  It  is  noted  again 
that  the  approach  taken  by  Hentschel  and  Procaccia  and  by  Kingdon  and  Ball  is  quite 
different  than  the  present  one.  It  is  especially  noteworthy  that  the  present  analysis  is  for  a 
free  shear  flow  rather  than  homogeneous,  isotopic  turbulence. 

Conclusions 

A  fractal  picture  of  constant  property  surfaces  in  round  jets  is  used  with  similarity 
arguments  and  empirical  data  to  obtain  an  expression  for  the  decay  of  mean  jet  fluid  mixture 
fraction  along  the  centerline  of  the  jet  which  is  of  appropriate  form  and  numerically  correct 
if  a  reasonable  value  for  Cd,  the  only  new  empirical  constant  in  the  modeling,  is  chosen. 

The  result  provides  support  for  the  hypothesis  that  constant  property  surfaces  in  free  shear 
flows  are  fractal  surfaces. 

For  Reynolds  number  independence,  the  R/  term  occurring  in  the  surface  area 
expression  must  cancel  the  R<  term  appearing  in  the  expression  for  ud  which  requires  that 
D  =  2  +  (p.  +  2)/6.  This  expression  for  D  is  equivalent  to  that  obtained  in  [4]  and  [5].  The 
observed  range  of  values  for  n  is  0.25  £  u  £  0.5  [4],  The  corresponding  range  for  D  is 
2.33  <  D  <  2.42,  a  result  which  compares  favorably  with  the  recent  experiments  of 
Sreenivasan  and  Meneveau  (31  for  free  shear  flows  and  with  the  results  of  Lovejoy  [21  for 


clouds. 


The  fractal  assumption  has  been  used  to  obtain  an  expression  for  mean  surface  area 

-  (10).  To  obtain  (10)  is  was  necessary  to  identify  an  outer  cutoff  area  and  to  assume  this 
area  is  equal  to  the  <Z>  surface  area.  A  well  defined  outer  cutoff  area  is  reasonable  for 
surfaces  which  fill  a  volume  that  is  relatively  narrow  in  one  direction  as  is  the  case  for  shear 
flows  but  is  most  likely  not  the  case  in  general.  The  assumption  relating  the  outer  cutoff 
area  to  the  <Z>  surface  is  justified  a  posteriori  by  the  success  of  the  modeling. 

Our  equation  for  D  can  be  obtained  without  identifying  an  expression  for  the  outer 
cutoff  area.  (15)  is  contingent  on  several  steps:  1)  expressing  dA^  in  terms  of  AJA0,  2) 

equating  A/Ag  to  (£j/£0)2'D  and  then  to  R/  l3/(4-u)]P  -  2](  3)  scaling  ud  by  u'R/1/2  -3/(4-u) 
with  u'  scaled  by  Ujd/x,  and  4)  requiring  Reynolds  number  independence.  1)  and  2)  form 
the  essential  features  of  the  fractal  hypothesis,  while  4)  is  supported  by  extensive 
experiment.  Of  the  four  steps  3)  is  perhaps  the  most  questionable.  As  discussed  above 
convective  fluxes  are  infinite  at  critical  points.  However  we  do  not  believe  these  points 
cause  problems  in  defining  average  fluxes.  As  an  alternate  to  the  Kolmogorov  scale  one 
could  choose  the  Taylor  microscale  to  represent  T)d.  An  analysis  similar  to  the  present  one  - 

-  without  correction  for  intermittency  -  which  uses  this  scale  gives  D  =  2  2/3.  This  result 
is  consistent  with  Mandelbrot's  for  Gauss- Kolmogorov  turbulence  [1]  but  is  inconsistent 
with  experiment  and  other  analyses. 

For  a  range  of  conditions,  chemical  reaction  in  turbulent  jet  diffusion  flames  occurs 
in  thin  sheets  referred  to  as  flamelets  [i»].  These  sheets  are  associated  with  constant 
mixture  fraction  surfaces,  i.e.,  the  stoichiometric  mixture  fraction  surface,  and  it  is  felt  that 
fractal  based  estimates  of  surface  area  will  be  very  useful  in  modeling  reacting  jets. 
Preliminary  work  in  this  direction  is  under  way  [12]. 

It  is  emphasized  that  the  fractal  picture  assumed  in  this  paper  is  open  to 
experimental  verification,  e.g.,  [3].  Laser  tomography  and  other  two-dimensional  imaging 
and  visualization  techniques  allow  one  to  study  the  fractal  character  of  curves  formed  by  the 


intersection  of  constant  property  surfaces  and  planes  of  illumination.  From  these  studies 
the  fractal  character  of  the  surfaces  themselves  can  be  inferred. 
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ABSTRACT 

A  technique  based  on  a  variation  of  Laser  Tomography  and  Laser 
Doppler  Velocimetry  has  been  developed  and  applied  to  measure  simultaneously 
the  flamelet  position  and  reactant  gas  velocity  in  pre  mixed,  turbulent  names.  The 
reactant  flow  is  seeded  with  a  silicon  oil  mist  which  rapidly  evaporates  and  burns 
in  the  flamelet.  Thus  the  intensity  of  scattered  light  from  a  laser  beam  decreases 
by  orders  of  magnitude  across  the  flamelet.  and  flamelet  position  along  a  line 
defined  by  the  beam  can  be  monitored  continuously.  The  signal  corresponding  to 
flamelet  position  is  differentiated  and  provides  information  on  flamelet  motion. 

Reported  here  are  pdfs  and  moments  of  flamelet  position  and  its 
derivative  for  methane-air  flames  at  two  different  equivalence  ratios.  These 
distributions  are  found  to  be  essentially  Gaussian,  and  it  is  observed  that  the  rms 
fluctuation  of  flamelet  position  increases  linearly  with  axial  distance.  Increasing 
equivalence  ratio  increases  these  rms  fluctuations,  but  spectral  analysis  shows  that 
the  spectral  characteristics  of  flamelet  motion  are  unaffected.  Further  comparison 
of  the  spectra  for  flamelet  motion  to  reactant  velocity  spectra  indicates  that  the 
flamelet  motion  is  controlled  by  the  turbulence  field  between  frequencies  of  300  to 
1500  Hr. 

Conditional  analysis  of  reactant  velocity  is  also  performed,  and  it  is 
found  that  rms  velocity  fluctuations  remain  nearly  constant  as  the  flamelet  is 
approached.  Higher  order  moments  of  the  conditioned  velocity  statistics  are  also 
reported. 

NOMENCLATURE 

A  Mean  reaction  rate  modelling  constant 

c  Mean  reaction  progress  variable 


<  w(  x)>  Mean  rate  of  chemical  reaction  per  unit  volume 
<*  Methane-Air  Equivalence  ratio 

a  Standard  deviation 

td  Detector  system  response 

I.  INTRODUCTION 

In  premixed  turbulent  flames,  for  mixtures  characterized  by  high 
activation  energies,  large  Damkohler  number,  and  for  moderate  levels  of 
turbulence,  the  major  heat  releasing  chemical  reactions  occur  in  thin  zones  which 
are  organized  into  sheets  [  1.2).  At  very  low  turbulence  levels,  these  sheets  can  be 
represented  by  unstrained  wrinkled-laminar-flames,  and  the  sheets  form  a  more  or 
less  contiguous  surface  in  the  flame  brush.  For  higher  levels  of  turbulence,  the 
structure  of  the  sheets  departs  from  that  of  an  unstrained  laminar  flame  as  the 
effects  of  sheet  curvature,  flame  stretch,  and  unsteady  flow  become  important. 
Even  so.  reaction  is  still  distributed  in  sheets,  and  the  dispersion  of  chemical 
reaction  into  small  volumes  of  highly  dispersed  fluid  occurs  only  at  high 
turbulence  levels.  These  sheets  of  chemical  reaction  are  referred  to  as  flamelets 
and  in  general  the  term  is  used  to  denote  both  the  high  temperature  reaction  zone, 
as  well  as  the  adjacent  convective-diffusive  layer  (the  preheat  zone). 

To  improve  our  understanding  of  turbulent  flames  and  to  support  model 
development,  the  geometry  of  flamelet  surfaces,  their  motion,  and  the  local 
structure  of  these  sheets  must  be  studied.  In  this  paper  we  report  on  a  technique 
for  simultaneous  measurements  of  flamelet  position  and  gas  velocity  developed 
primarily  to  measure  gas  velocities  conditioned  on  distance  from  the  flamelet 
surface.  The  approach  also  allows  one  to  obtain  information  on  the  statistics  of 
flamelet  position  in  space  and  flamelet  motion.  These  measurements,  for  two 
different  flames,  are  also  reported  and  discussed. 


Instantaneous  droplet  diameter 
Initial  droplet  diameter 
Burning  constant  for  oil  droplets 
Turbulent  integral  length  scale 
Oil  mist  droplet  concentration 
Laser  Power 

Flamelet  position  probability 

Turbulent  Reynolds  number 

Space  coordinate  along  laser  beam 

Measured  photomultiplier  tube  signal 

Smamwise  velocity  component 

Cron  sauam  velocity  component 

EnaemMe  average  volumetric  chemical  reaction  rate 


Still  photography,  high-speed  cinematography,  shadow  and  Schlieren 
techniques  have  all  been  used  to  obtain  information  on  flamelet  surfaces.  These 
techniques  have  proven  useful,  but  their  interpretation  is  ambiguous  primarily 
because  they  are  all  line-of-sight  techniques.  Hence  they  provide  path  integrated 
information. 

Laser  tomography  and  related  flow  visualization  techniques  avoid  line-of- 
sight  problems  and  are  potentially  much  more  useful  for  studying  flamelet  surface 
characteristics  than  the  above  methods.  Both  still  photographs  and  cine-films  can 
be  made,  thus  allowing  flamelet  dynamics  to  be  investigated.  However,  it  should 
be  noted  that  all  these  methods  are  subject  to  the  inherent  limitation  of  visualizing 
a  curve  generated  by  the  intersection  of  the  flamelet  surface  and  a  plane,  rather 
than  the  flamelet  surface  itself. 

In  laser  tomography,  as  applied  to  premised  flames,  the  reactant  flow  is 
seeded  with  an  oil  mist  and  the  flow  is  illuminated  with  a  sheet  of  laser  light. 
Since  the  oil  rapidly  vaporizes  in  the  reaction  zone  of  the  flamelet.  scattering  from 
the  laser  sheet  is  observed  only  from  the  reactants,  and  the  flamelet  is  marked  by 
the  boundary  between  the  scattering  and  non-scattering  regions.  This  boundary  is 
easily  recorded  photographically.  Several  investigators  have  used  tomography 
for  visualization,  including  its  originator,  Boyer  (3),  and  Borghl  (4). 

A  variation  of  laser  tomography  (5]  allows  for  the  measurement  of 
flamelet  location  along  a  line  in  space,  defined  by  the  laser  sheet  illumination  and 
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(he  image  of  a  slit  on  this  sheet.  In  this  technique,  a  photomultiplier  tube 
measures  the  intensity  of  scattered  radiation  from  the  line  in  space,  where  the  slit 
is  placed  in  front  of  the  photomultiplier  tube,  and  scattered  photons  are  collected 
by  a  lens  which  in  mm  makes  an  image  of  die  laser  illumination  shot  on  this  slit 

The  integrated  intensity  of  scattered  light  falling  on  the  detector  surface 
depends  on  the  position  of  the  flamelet  along  the  line,  and  thus  the  position  can  be 
monitored  by  recording  the  photomultiplier  output 

Boyer  et  al.  (5|  used  this  position  measurement  technique  along  with  laser 
doppler  velocimetry  (LDV)  to  study  flamelet  motions  and  gas  velocities  in  very 
weak  turbulence.  They  record  for  analysis  three  signals:  the  flamelet  position 
signal  and  its  derivative,  and  the  LDV  signal  from  the  measurement  of  a  single 
velocity  component  at  a  fixed  point  in  space.  The  data  are  analyzed  for  the  power 
spectra  of  each  signal  in  a  frequency  range  up  to  10  Hz.  A  single  flame  condition 
is  studied  and  only  spectral  data  are  reported  and  discussed.  The  derivative  of  the 
flamelet  position  signal  is  treated  as  a  component  of  the  flamelet  velocity.  As 
discussed  below,  this  treatment  of  the  derivative  signal  is  subject  to  interpretation. 

Suzuki  and  Mirano  (6|  have  reported  measurements  of  flamelet  orientation 
and  speed  in  a  plane,  obtained  at  a  point  using  a  triple  sensor  ion  probe.  Their 
analysis  of  data  assumes  that  the  flamelet  may  be  represented  by  a  plane  which  is 
parallel  to  the  axis  of  their  probe,  and  that  the  spatial  orientation  of  the  plane 
remains  constant  as  it  crosses  the  probe  volume.  It  is  claimed  that  highly  curved 
flamelets.  and  flamelets  which  cross  at  a  large  angle  to  the  probe  axis  have 
characteristic  ion  current  signatures  which  allow  for  these  events  to  be  recognized 
and  discriminated  against.  At  best  this  procedure  will  bias  the  results.  Very  high 
instantaneous  flamelet  surface  speeds  are  reported,  which  raises  questions  about 
the  assumptions  made  in  the  data  analysis  and  about  Suzuki  and  Hirano's  ability 
to  recognize  flamelet  crossing  events  which  violate  their  assumptions. 

To  date  the  above  work  and  related  references  are  the  only  reports  in  the 
literature  where  quantitative  measurements  of  flamelet  position,  orientation  and 
motion  are  reported.  Boyer  et  al.  (3)  record  data  which  could  be  used  to  make 
conditioned  velocity  measurements,  but  no  results  of  this  type  are  reported. 

In  this  paper  we  describe  a  technique  for  making  flamelet  position  and  gas 
velocity  measurement  which  is  similar  to  that  of  Boyer  et  al.  [3  j.  The  method  is 
applied  to  higher  velocity,  more  turbulent  flows  than  that  of  Boyer  et  al.  and  the 
objectives  are  quite  different  The  laser  tomographic  technique  is  used  to  measure 
flamelet  position  along  a  laser  beam;  gas  velocity  at  a  fixed  point  on  the  laser 
beam  is  obtained  by  LDV.  The  flamelet  position  signal,  its  derivative,  and  the  gas 
velocity  signal  are  recorded.  These  data  records  are  then  analyzed  for  velocity 
statistics  conditioned  on  the  distance  of  the  velocity  measurement  point  from  the 
flamelet  It  is  expected  that  effects  of  combustion  on  the  turbulence  will  be 
evident  in  such  data.  The  statistics  of  the  flamelet  position  along  the  laser  beam 
are  also  of  interest  and  have  implications  for  the  nature  of  the  flamelet  motion. 
For  instance,  the  distribution  will  be  bi-modal  if  large  amplitude  sine  waves  are 
causing  flamelet  motion.  Finally,  time  derivatives  of  flamelet  position  are 
measured  and  compared  with  gas  velocity  data. 


II.  EXPERIMENT 
A.  Apparatus 

1.  Burner  Facility.  Figure  1  is  a  schematic  of  the  burner  exit  and 
coordinate  system.  The  turbulent  flame  is  stabilized  by  a  1.3  mm  diameter  rod 
mounted  across  the  exit  plane  of  the  30  mm  diameter  cylindrical  burner. 
Turbulence  is  generated  by  interchangeable  woven  wire  screens  (10  mesh  and  8 
mesh)  which  are  mounted  in  the  central  air-fuel  jet.  1.0  cm  below  the  exit  plane. 
The  fuel  is  commercial  grade  methane  supplied  from  high  pressure  bottles.  Air. 
filtered  and  dried,  is  supplied  from  high  pressure  tanks  which  are  charged  by 
compressor  (both  gases  are  metered  by  needle  valves,  and  flow  rates  are 
monitored  with  rotameters).  The  burner  is  equipped  with  an  annular  jet  whose 
velocity  may  be  matched  to  that  of  the  inner  jet  in  order  to  shield  the  flame  from 
the-mixing  layer  formed  with  room  air.  Due  to  the  limited  air  Stonge  capabiltiy, 
and  a  primary  desire  to  use  long  averaging  times  in  order  to  more  closely 
approximate  the  higher  order  momenta  of  the  measured  quantities,  this  feature 
was  not  utilized.  All  measurements,  however,  are  perforated  well  away  from  the 
edge  of  the  jet.  and  it  is  felt  that  the  effect  of  the  mixing  layer  in  this  region  is 
negligible.  This  burner  has  been  used  in  several  previous  studies  at  Coraeil,  e  g. 
(7)  and  (8). 

2.  LDV  System.  Reactant  gas  velocities  tre  measured  using  in  LDV 
apparatus  consisting  of  an  argon-ion  laatr  (Lexet,  Modal  93),  commercial  optics 
with  a  Bragg  cell  frequency  shifter  (Thermosystams,  Inc.,  900  series)  and  a 
counter  type  signal  processor  (Thermotystems.  !nc.  Model  1980).  A  dual  beam, 
teal  fringe  optical  configuration  is  employed.  The  plane  containing  the  two  beams 
is  aligned  such  that  the  major  aais  of  the  ellipsoidal  measurement  volume  lies 
parallel  to  the  flame  stabilizer  rod.  Scattered  light  from  the  seed  particles  is 
collected  at  a  90*  angle  from  the  incident  radiation,  and  a  0.23  mm  diameter 
pinhole  at  the  entrance  to  the  photomultiplier  tube  limits  the  major  axis  diameter  to 
about  0.8  nun.  white  the  minor  axis  diameters  are  approximately  0.13  mm. 

Several  checks  on  the  validity  of  the  velocity  measure  menu  were 
performed  on  both  a  laminar  and  turbulent  jet  in  an  attempt  to  maximize  the  data 
rate  without  significantly  distorting  the  messured  velocity  probability  density 
function  (pdf).  Laminar  jet  measurements  ware  made  mainly  to  verify  the  correct 
operation  of  the  Bragg  cell.  In  two  independent  data  seta,  the  measured  mean 
latarai  velocity  was  0JJ  m/sec.  but  one  of  the  data  seta  showed  secondary  peaks  in 


Figure  1.  Schematic  of  Burner  and  Coordinate  System. 

the  tails  of  the  pdfs,  which  were  present  at  all  data  rates.  This  does  not  appear  in 
the  streamwise  velocity  component  pdfs,  nor  in  all  of  the  lateral  component 
measurements,  and  has  been  attributed  to  intermittent  leaking'  of  secondary 
frequencies  in  the  downmix  circuit  of  the  Bragg  cell.  Tests  performed  on 
turbulent  jets  demonstrated  that  data  rates  of  about  10  kHz  were  feasible  without 
significantly  broadening  or  distorting  streamwise  component  velocity  pdfs. 
Above  10  kHz  the  pdfs  begin  to  broaden  and  depan  significantly  from  a 
Gaussian.  Cross-stream  component  pdfs  were  qualitatively  similar,  again 
showing  little  broadening  until  dau  rates  were  pushed  above  10  kHz,  although  an 
occasional  disconcerting  shift  in  the  measured  mean  velocity  was  observed.  The 
broadening  of  the  pdf  as  photomultiplier  tube  pain  and/or  amplifier  gain  is 
increased  is  to  be  expected:  with  increased  gain,  marginal  particles  passing 
through  the  edges  of  the  measurement  volume  are  seen.  In  these  regions  the 
fringes  are  not  as  parallel  and  planar  as  they  are  near  the  center  of  the 
measurement  volume,  resulting  in  an  increase  in  the  measured  variance.  The 
occasional  shift  in  measured  mein  velocity  is  more  difficult  to  explain;  again 
suspicion  falls  on  the  downmix  circuit.  Velocity  data  was  taken  at  data  rates 
corresponding  to  10  kHz  when  the  measurement  volume  is  wholly  within  the 
reactants.  The  preliminary  results  presented  in  this  paper  are  repotted  with  the 
understanding  that,  though  qualitatively  correct,  more  checks  of  cross-stream 
velocity  measurement  validity  are  required,  and  are  being  performed  in  the  on¬ 
going  work. 


Figure  2.  Schematic  of  Experiments!  and  Data  Acquisition  System. 


3.  Laser  Tomographic  System.  A  sketch  of  the  liser  tomographic 
system,  as  well  as  the  relative  position  of  the  LDV  system,  is  shown  in  figure  2. 
A  6  mW,  He-Ne  laser  beam  la  pasted  over  the  center  of  the  burner  exit  normal  to 
the  vertical  plane  containing  the  flame  stabilizer  rod.  A  plano-convex  lens  (f  * 
1000  m,  f  «  30  mm)  is  used  to  reduce  the  l/*1  diameter  of  the  beam  to  about  0.53 


mm  where  it  passes  over  the  burner  exit.  A  poruon  of  this  beam  is  imaged  by  a 
symmetric,  biconvex  lens  (f  -  127  mm,  *  -  50  mm)  onto  the  entrance  sPit  of  the 
photomultiplier  tube  with  a  magnification  equal  to  1.0.  The  entrance  slit  has  a 
vertical  dimension  of  0J  mm  and  a  length  of  9.66  mm;  thus,  accounting  for  the 
angle  between  the  axis  of  the  collection  optics  and  the  laser  beam,  1.0  cm  of  the 
beam  is  imaged  onto  the  entrance  slit. 

The  signal  out  of  the  photomultiplier  tube  is  amplified,  low-pass  filtered  at 
J  kHz.  and  differentiated.  The  analog  differentiating  circuit  |9)  was  constructed 
with  a  high  slew  rate  op-amp  ind  precision  components.  The  component  values 
result  in  cut-off  frequencies  of  26  kHz.  to  which  point  the  circuit  differentiates, 
and  33  kHz.  after  which  the  circuit  integrates  to  avoid  'ringing'  caused  by  high 
frequency  signal  components.  The  circuit  was  calibrated,  and  found  to  follow  the 
theoretical,  first  order  gain  curve  very  closely  up  to  about  3  kHz.  where  the  gain 
begins  to  roll  off  gradually.  The  response  of  the  circuit  to  a  pure  sine  wave  at  3 
kHz  was  about  8.3%  below  the  calculated  time  derivative.  Data  analysis  was 
performed  using  the  theoretical,  linear  gain  curve. 

Seeding  of  the  reactants  for  both  laser  tomography  and  LDV 
measurements  is  accomplished  using  atomizing  nebulizers  similar  to  those  used 
by  Cheng  (101  *"d  described  by  Durst  Melling  and  Whitelaw  {ill,  who  state  that 
with  careful  design  a  mean  seed  particle  diameter  of  near  1.0  micron  can  be 
achieved.  Relatively  clean  Doppler  bursts  from  the  LDV  scattering  volume 
indicate  that  the  majority  of  seed  panicles  are  less  than  about  3  microns  in 
diameter,  though,  as  noted  below,  particles  as  large  as  100  microns  are  present  as 
well.  Dow  Coming  50  centistoke  silicon  oil  is  used.  The  seed  is  introduced  in 
the  mixing  chamber  in  the  base  of  the  burner,  where  it  is  thoroughly  mixed  with 
the  reactants.  The  mixture  then  passes  through  several  fine  wire  mesh  screens, 
which  serve  to  dampen  any  large  scale  velocity  fluctuations  in  the  flow  and 
disperse  the  seed  uniformly  before  it  passes  to  the  exit  section  of  the  burner.  The 
particle  number  density  at  the  exit  plane  of  the  burner  is  estimated  to  be  about 
1000  particles  per  mm}. 

All  three  signals;  gas  velocity,  flamelet  position  and  its  derivative  are 
convened  to  digital  format  (Date!,  ST-PDP  Series)  and  stored  for  analysis  on  a 
DEC  PDP  11-34  computer. 

D.  Details  of  Flamelet  Position  Measurement 

The  integrated  intensity  of  scattered  light  from  a  centimeter  of  the  laser 
beam  is  measured  by  a  photomultiplier  tube.  As  outlined  above,  scattered  light  is 
observed  only  from  the  reactants,  and.  subject  to  the  following  conditions,  the 
photomultiplier  tube  output  is  proposed  to  be  linearly  related  to  the  position  of  the 
flamelet  along  the  laser  beam;  (1)  the  seed  distribution  in  the  reactants  is  uniform, 

(2)  overall  scattering  is  not  sufficient  to  significantly  auenuate  die  laser  intensity, 

(3)  the  detector  sensitivity  is  uniform  along  the  beam,  (4)  the  beam  has  a 
negligibly  small  diameter,  and  (5)  the  evaporation/  combustion  time  of  the  oil 
droplets  is  small  as  compared  to  their  passage  time  through  the  flamelet.  These 
conditions  will  now  be  discussed. 

The  photomultiplier  output  may  be  written  as  the  integral  of  the  product  of 
the  laser  beam  intensity,  the  oil  mist  concentration  (n).  an  average  scattering 
cross-section  tor  the  mist  (o),  and  the  local  detector  system  response  (td),  which 
in  general  is  a  function  of  position.  The  limits  of  integration  are  the  volume 
formed  by  the  laser  beam  and  the  slit  image  on  the  laser  beam.  Since  scattering  in 
the  product'  is  negligibly  small,  integration  may  be  limited  to  those  portions  of 
the  mea.iur  <tent  volume  containing  reactants.  The  integral  is  evaluated  first  by 
integrating  over  the  laser  beam  cross-section  leading  to  an  expression  in  terms  of 
quantities  averaged  over  the  cross-section; 

S  -  Jnax^Pds. 

where  S  is  the  measured  signal,  s  is  the  distance  along  the  laser  beam,  and  P  is 
the  laser  power.  Integration  is  along  portions  of  the  beam  containing  reactants, 
for  uniform  seed  distrioution  in  the  reactants,  no  is  constant;  tram  above,  %  and 
P  are  assumed  constant.  Therefore  S  is  proportional  to  the  length  of  the  laser 
beam  subtending  the  reactants,  and  if  the  flamelet  cresses  the  beam  only  once  in 
the  field  of  view,  S  is  proportional  to  the  position  of  the  flamelet  along  the  beam. 

So  far.  the  angle  of  intersection  between  the  laser  beam  and  the  flamelet 
has  been  ignored.  For  a  planar  flamelet  and  a  cylindrical  laser  beam  it  is  easily 
seen  that  variations  in  the  angle  of  intersection  do  not  influence  the  measurement 
unul  the  beam  and  surface  approach  •  parallel  state.  Except  for  this  extreme,  the 
measurement  gives  the  point  of  intersection  of  the  flamelet  and  the  axis  of  the 
laser  beam.  Flamelet  curvature  will  affect  this  relationship  and  therefore  it  is  best 
to  keep  the  beam  diameter  to  e  minimum,  thereby  minimizing  curvature  effects. 

The  position  measurement  system  is  calibrated  by  removing  the  turbulence 
generating  grid  and  holding  a  laminar  flame,  *  -  0.8.  on  the  burner.  The  flame  is 
then  traversed  along  the  laser  beam,  while  recording  the  output  of  the 
photomultiplier  tube  as  a  function  of  flame  position.  A  least  squares  linear  curve 
fit  is  applied,  resulting  in  an  expression  for  flamelet  position  as  a  function  of 
photomultiplier  tube  output  This  expression  is  accurate  to  within  ±0.3  mm,  and 
to  within  approximately  ±0.13  mm  within  the  central  8  mm  of  the  measurement 
volume.  System  response  was  found  to  be  quite  linear,  substantiating  the 
assumptions  ot  umtorm  no  and  constant  %  (t-igure  i). 

An  estimate  of  droplet  aviporxoon/combusbon  rime  of  approximately  10  * 
seconds  ha*  been  reported  by  8oyer  (3);  our  malysrs  using  the  so-called  *d‘  lew’ 
of  droplet  combustion  (11)  confirm*  this  report  for  droplets  up  to  approximately 


Figure  3.  Position  Signal  Calibration  Curve. 

10  pm  in  diameter.  For  a  typical  velocity  component  normal  to  the  flame  front  of 
0.30  m/s,  and  a  flamelet  thickness  of  about  1.0  mm,  the  residence  time  of  the 
droplet  within  the  flamelet  is  about  3.3  x  10-J  seconds.  Experimental  data, 
however,  obtained  from  point  measurements  of  scattered  light  (spatial  resolution 
of  about  0.15  mm)  near  the  flamelet.  indicate  that  some  droplets  persist  after 
passing  through  the  flamelet.  Only  after  a  distance  of  about  1.5  mm  from  the 
flamelet  have  they  completely  disappeared.  This  would  indicate  thet  some  very 
large  droplets  are  present  in  the  seeded  reactants,  perhaps  as  large  is  1 00  pm  in 
diameter.  These  particles  are  rare,  and.  due  to  the  calibration  procedure  described 
above,  their  existence  in  no  way  biases  our  measurements. 

An  additional  concern  related  to  droplet  evaporation  time  is  die  effect  of 
increased  equivalence  ratio,  resulting  in  increased  flame  temperature,  on  the 
flamelet  position  measurements.  Since  the  burning  constant,  K,  in  the  d1  law.  d1 
•  do*  -  Kt,  depends  logarithmically  on  the  flame  temperature,  this  ettect  is  slight 
Our  estimates  indicate  that  the  change  in  measured  flamelet  position  due  to 
increased  equivalence  ratio  (♦  «  0.8  to  ♦  -  1.0)  is  about  lO^mm.  Finally,  the 
heat  release  due  lo  oil  droplet  combustion,  less  the  energy  required  to  completely 
vaporize  the  droplet  has  been  compared  to  the  heat  release  of  the  methane/air 
combustion.  A  ratio  of  I  to  200  indicates  (hat  this  effect  is  negligible. 

Spatial  resolution  of  the  position  measurement  is  limited  by  noise  in  the 
signal.  Possible  significant  sources  of  noise  are  photon  shot  noise,  large  scale, 
low  frequency  fluctuations  in  seed  concentration,  and  marker  shot  noise  (i.e., 
random  fluctuations  in  seed  concentration  (13);  the  latter  being  by  far  the  most 
significant  problem.  Scattering  signals  from  single  particles  ire  quite  high,  and 
photon  shot  noise  is  not  a  problem.  Two  point  time  correlation  measurements 
with  the  points  lying  on  a  streamline  show  zero  correlation  between  signals  from 
pairs  of  points  except  for  a  very  small  value  at  a  lag  time  corresponding  to  the 
flow  time  between  the  points  (Figure  4).  Thus  large  scale  fluctuations  in  seed 
flow  are  negligible  compared  to  other  noise  sources,  namely  marker  shot  noise. 

For  identical  particles,  marker  shot  noise  would  be  expected  to  follow 
Poisson  statistics.  Variations  in  particle  scattering  cross-section  at  low  particle 
arrival  rates  might  cause  non-Poisson  behavior.  However,  measurements  of 
noise  pdfs  show  Ctussisn  distributions  as  would  be  expected  for  Poisson 
processes  and  large  samples.  Marker  shot  noise  can  be  reduced  by  increasing  the 
number  of  particles  in  the  measurement  volume  and  by  reducing  the  band  width 
of  rite  detector. 

To  increase  particle  loadings,  four  atomizing  nebulizers  are  connected  in 
parallel,  and  a  large  fraction  of  the  air  flow  passes  through  them.  Higher  particle 
concentration  levels  are  possible  with  condensation  type  seeders  (which  give 
smaller  particles  as  well);  this  type  of  device  should  be  explored.  An  increase  in 
measurement  volume  will  reduce  marker  shot  noise  with  a  concomitant  loss  in 
spatial  resolution.  Limiting  band  width  reduces  temporal  resolution  as  well  as 
noise.  Taylor's  hypothesis  may  be  used  to  express  temporal  resolution  in  spatial 
terms.  For  U  -  3.8  m/sec  and  a  temporal  resolution  of  0.1  ms  (the  minimum 
time  between  successive  sample*  of  the  data  channels),  the  equivalent  spatial 
resolution  is  0.38  mm,  which  is  of  the  same  order  as  the  nominal  diameter  of  the 
laser  beam.  Thus  the  two  are  consistent  with  each  other.  In  the  measurements 
reported,  the  position  signal  is  low  pats  filtered  at  5  kHz  before  storage.  With 
this  filter  setting,  a  1.0  cm  long  scattering  volume,  and  the  0.35  mm  laser  beam 
diameter,  the  measured  signal-to-noiae  ratio  is  43  at  full  intensity  (the 
measurement  volume  being  entirely  in  the  reactants)  and  37  at  half  intensity.  For 
a  Gaussian  distribution  of  noil*,  theta  notae  figures  correspond  to  ipatial 
uncertainties  at  a  95%  confidence  level  of  ±0.4  mm  and  ±03  mm  for  full  Kale 


« 


and  half  scale  measurements,  respectively.  These  numbers  compare  favorably 
with  the  spatial  and  temporal  resolutions  of  the  measurement  system,  which 
confirms  our  choice  of  these  two  quantities. 


time  (tee) 

Figure  4.  Two-Point.  Time  Conelatioa  of  Scattering  Signal. 


C.  Data  Analysis 

The  three  dare  signals — flamelet  position,  its  derivative,  and  the  reactant 
velocity —  are  analyzed  in  several  different  ways:  for  pdfs,  power  spectra,  and 
certain  conditioned  satisfies.  Standard  computer  subroutines  are  used  for  most  of 
these  analyses.  However,  some  comments  are  in  order  regarding  dau  analysts 
and  interpretation. 

It  is  tempting  to  interpret  the  derivative  of  the  position  signal  as  the 
flamelet  velocity,  and  for  those  cases  when  the  flamelet  is  perpendicular  to  the 
laser  beam  this  is  a  valid,  unambiguous  interpretation.  In  general,  however,  this 
view  must  be  scrutinized.  It  is  quite  clear  that  the  derivative,  in  general,  does  not 
give  the  velocity  of  a  fixed  point  on  the  flamelet  while  it  does  give  the  rate  at 
which  the  flamelet  position  along  a  line  changes  in  time.  Only  velocity 
componena  normal  to  the  flamelet  surface  will  cause  it  to  move.  Tangential 
velocity  componena  move  fluid  along  the  surface  but  do  not  move  the  surface 
itself,  while  gradiena  in  the  normal  velocity  cause  wrinkling  and  rotation  of  the 
surface.  Thus  surface  points  where  the  surface  lies  perpendicular  to  the  laser 
beam  move  parallel  to  the  beam.  Otherwise,  they  move  at  an  angle  to  the  beam 
and  the  motion  of  the  surface  along  the  laser  beam  does  not  correspond  to  the 
motion  of  a  fixed  surface  point  In  general,  when  one  refers  to  the  velocity  of  a 
surface,  it  is  in  regard  to  a  fixed  point  on  the  surface.  We  do  not  measure  a 
component  of  this  velocity  here.  However,  the  derivative  quantity  we  do  measure 
does  given  information  on  the  flamelet  motion  in  space  and  therefore  is  of 
interest 

A  major  objective  of  this  work  is  to  detect  effects  of  the  flamelet  on  the 
turbulence  field.  This  may  be  done  in  part  by  conditioned  velocity  measurements. 
In  the  past  velocity  measurements  conditioned  on  being  either  in  reactants  or 
producu  have  been  reported  (14,15).  As  is  seen  below,  such  measure  menu  do 
not  fully  reveal  the  turbulence  effects.  What  is  needed  is  a  conditioning  of  the 
velocity  by  distance  away  from  the  instantaneous  flamelet  position.  The 
measurements  reported  here  are  a  first  step  in  that  direction. 

Using  the  flamelet  position  signal,  velocity  data  are  segregated  by  the 
distance  along  the  laser  beam  from  the  flamelet  to  the  measurement  point  For 
this  segregation,  the  distance  is  divided  into  I  mm  segments,  i.e..  0  to  1  mm 
away.  I  to  2  mm  away,  etc.  Once  segregated,  the  velocity  data  are  used  to 
determine  means,  variance,  skewness  and  lturtosis  of  the  conditioned  velocity 
distributions. 


III.  RESULTS  AND  DISCUSSION 

Dau  was  taken  for  methane-air  flames  at  equivalence  ratios  of  0  J  and 
1.0.  For  each  flame,  measurements  were  performed  at  four  different  axial 
locations:  3J.  6.0,  6.3  and  7.0  cm  above  the  turbulence  generating  grid.  Two 
lea  of  daid  were  taken  at  eaeh  location,  one  'slow'  set  at  a  sampling  rat*  equal  to 
twice  the  integral  time  scale  (based  on  Taylor's  hypothesis  and  measured  integral 
length  scales)  from  which  averages  and  moms  ms  as  computed,  and  one  ’fast’  set 
which  ts  used  for  spectral  analysis.  Each  slow  dau  set  consists  of  40.960  points. 
These  points,  by  virtue  of  the  slow  sampling  rate.  «e  independent  (161  and. 
assuming  s  Gaussian  pdf.  result  in  convtrganca  of  the  calculated  variance  to 
within  0.7%  of  dM  true  me  sewed  variance.  Skewness  and  kernels  converge  n 


within  1.3  and  2.9%,  respectively.  The  fast  data  sets  consist  of  102.400  points 
sampled  at  10  kHz,  in  blocks  of  duration  0. 1  sec.  As  the  flamelet  position  signal 
is  filtered  at  5  kHz.  this  should  result  in  alias-free  power  spectra. 

Table  I  shows  the  turbulence  properties  and  scaling  of  the  approach  flow 
as  determined  by  previous  measure  menu  using  hot-wire  anemometry  [17],  and  is 
presented  here  for  reference. 

In  Table  2  are  shown  the  Tint  four  moments  of  the  flamelet  position 
signal.  These  clearly  indicate  both  the  effect  of  equivalence  ratio  on  the 
magnitude  of  the  flamelet  fluctuations  and  the  Gaussian  nature  of  these 
fluctuations.  The  increase  in  the  rmt  fluctuations,  which  may  be  taken  as  a 
measure  of  the  turbulent  flame  brush  thickness,  a  the  equivalence  ratio  is 
increased  is  similar  to  the  results  recently  reported  by  Shephard  ct  al.  [181,  who 
have  attributed  it  to  increased  heat  release  in  the  reaction  zone.  The  very  nearly 
linear  increase  of  these  fluctuations  with  axial  distance  should  also  be  noted.  The 
Gaussian  nature  again  is  similar  to  recently  repotted  resuiu  (1 8. 1 9],  and  further 
substantiates  the  assumption  of  a  Gaussian  pdf  in  the  model  developed  by 
Namazian,  et  al.  (20)  to  determine  mean  density  profiles  across  the  flame  brush. 
Figure  5  shows  a  representative  of  the  measured  pdfs. 

Table  l 

Turbulence  Fropertite  and  Seeling 


U  -  3.8  m/sec 

8  mesh  arid 

x(cm) 

5.5 

6.0 

6.5  7.0 

u/U  (%) 

3.5 

3.3 

3.2  3.1 

l(mm) 

2.48 

2.53 

2.64  2.70 

Rt(ut/v) 

24.4 

23.5 

23.9  23.5 

Table  2 

Moments 

of  Flamelet  Position  Signal 

y  (mm) 

rms  (mm)  skew  kurtosis 

♦  nfl.8: 
x  -  5.5  cm 

4.02 

0.77 

-0.08 

2.82 

x  -  6.0  cm 

4.40 

0.85 

-0.09 

2.84 

x  -  6.5  cm 

4.95 

0.93 

-0.12 

2.80 

x  »  7.0  cm 

5.49 

1.00 

-0.12 

2.81 

6  m  10: 

X  -  5.5  cm 

6.55 

0.74 

0.05 

3.07 

x  -  6.0  cm 

7.85 

0.95 

0.10 

2.83 

x  -  6.5  cm 

8.87 

1.06 

0.08 

2.87 

x«  7.0  cm 

10.50 

1.39 

0.07 

2.68 

Figure  5.  Flamelet  Position  PDFs  for  6  -  0.8. 

Axial  Positions  of  5 S,  6.0, 6.5  and  7.0  cm. 


It  Is  tlso  interesting  lo  consider  the  relationship  of  the  pdf  of  flamelet 
position.  P(x),  to  the  chemical  closure  problem.  If  <«o(x)>  is  the  mean  rate  of 
chemical  reaction  per  unit  volume  at  a  given  x.  and  atf  is  the  ensemble  average 
volumetric  chemical  reaction  raw  per  flamelet  crossing,  then 

<««*)»  -  o^n*)  (•) 


1 


Figure  6.  Comparison  of  P(x)  and  Ac(l-c). 


Frequently  (Bray  et  al.  (21)),  <o»  is  modelled  as  a  function  of  c(a).  the 
mean  reaction  progress  variable:  «u>  a  c(l-c).  It  is  easily  shown  that  for  thin 
flame  lets. 


era)  -  J  P(x)  dx'  (b) 


Using  this  relationship,  one  can  compare  the  spatial  variation  of  <o» 
predicted  by  c(I-c)  and  that  given  by  Pfx).  To  do  so,  P(x),  c(x)  and  c(l-c)  can 
be  obtained  as  a  function  of  x  using  the  above  relationships  and  the  measured 
Pfx).  We  have  made  such  a  comparison  using  a  normal  distribution  to  represent 
the  measured  Pfx).  primarily  for  convenience.  The  results  are  shown  in  Figure  6, 
where  P  and  Ac(l-c)  are  plotted  against  c.  A  is  chosen  to  match  Acfl  -c)  and 
Pfx)  at  c  -  05;  A  -  1.596.  It  can  be  seen  that  Acfl-c)  gives  an  extremely  good 
approximation  of  P  in  c  space,  which  accounts  for  its  success  in  modeling  mean 
combustion  rates.  The  points  labeled  o.  2 a  and  3o  in  the  figure  indicate  c  values 
associated  with  points  1, 2  and  3o  frms)  away  from  the  center  of  the  flame  brush 
in  physical  space. 

Table  3  shows  the  mo  menu  of  the  derivative  of  the  flamelet  position 
signal.  Notable  features  are  the  rms  values  again  increasing  with  both  axial 
distance  and  equivalence  ratio. 

The  spectra  of  the  flamelet  position  signal  and  the  spectra  of  its  derivative 
are  plotted  in  Figures  7  and  8.  The  flamelet  position  spectra  for  ♦  -  0.8  »nd  ♦  - 
1 .0  are  very  similar,  adding  further  support  to  the  assertion  ( 1 8]  that  for  the  same 
inlet  turbulence  conditions  an  increase  in  equivalence  ratio  increases  only  the 
amplitude  of  the  flamelet  motion,  while  the  spectral  characteristics  we  unchanged. 
Comparison  of  the  spectra  for  the  time  derivative  of  flamelet  position  shows  that 
they.  too.  remain  quite  similar  as  the  equivalence  ratio  increases.  This  is  to  be 
expected  since  these  spectra  can  be  derived  from  the  corresponding  flamelet 
position  spectra.  A  comment  on  the  shape  of  these  spectra  is  in  order.  It  is  well 


Table  3 

Moments  of  the  Derivative  of  the 


Flamelet  Position 

Signal 

rms  (m/sac) 

skew 

kunotu 

*  -  0.8: 
x  -  5  5  cm 

2.44 

0.21 

357 

x  -  6.0  cm 

2.54 

0.29 

4.11 

x  -  6.5  cm 

2.63 

0.32 

4.37 

x  -7.0  cm 

2.67 

0.37 

4.7* 

4-  1.0: 
x  •  5-5  cm 

2.78 

•0.36 

J23 

x  -  6.0  cm 

3.4J 

•0.50 

5.38 

x  -  6.5  cm 

3.51 

-0.77 

6.17 

x  •  7.0  cm 

4.31 

-0.91 

640 

for  p  -  0.8  and  ♦  -  1.0. 


Figure  8.  Comparison  of  Spectra  for  the  Derivative 
of  Flamelet  Position;  ♦  -  0.8  snd  ♦  -  1.0. 


known  [16)  that  the  value  of  a  one-dimensional  frequency  spectra  of  some 
physical  process  at  zero  frequency  is  propoitiontl  to  the  integral  lime  scale  of  that 
process.  Since  the  fltmelet  position  is  a  statistically  stationary  process,  the 
integral  scale  of  its  derivsuve  must  vanish.  This  is  borne  out  by  the  shape  of  the 
spectre  for  the  derivative  of  flamelet  position. 

In  Figure  9  the  spectrum  for  the  derivative  of  flamelet  position  is 
compared  with  the  one-dimensional  v -component  velocity  spectra.  The  velocity 
spectre  is  obtained  by  LDV  measurements  at  a  mean  data  raw  of  10  kHx.  snd  can 
therefore  be  expected  to  be  accurate  up  to  frequencies  of  2  or  3  kHt.  at  which 
point  the  discrete  nature  of  the  counter's  snalog  output  signal  causes  distortion. 
As  can  be  seen,  the  spectre  an  nearly  identical  between  frequencies  of  about  500 
to  1500  Hx.  This  type  of  tom  parson  has  siso  been  reported  by  Shepherd  et  al. 
(18),  though  a  coincidental  region  between  100  and  1000  Hx  was  reported  The 
discrepancy  can  probably  be  txplained  by  the  different  inlet  turbulence 
characteristics,  snd  it  appears  reasonable  to  assert  that  within  this  spectral  region 
the  fltmelet  movement  is  controlled  by  the  turbulence  field.  Moreover,  at  low 
frequencies  the  velocity  spectre  is  sffected  by  aliasing;  the  energy  of  high 
frequency  components  whose  velocity  vector  is  not  aligned  with  the  direction  of 
measniement  is  alitsod  to  low  frequencies.  Thus,  although  the  spectra  diverge  at 
low  frequencies,  h  is  quit*  plausible  that  flamelet  morion  in  this  spectral  region  is 
siso  turbulence  controlled. 


V. 

VJ 

v' 
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derivative  o<  flamelet  petition 


Frequency  (Hz) 

Figure  9_  Comparison  of  Spectra  for  the  Derivative  of  Flamelet 
Position  and  One -dimensional  V -component 
Velocity  Spectra.  9  •  0.8. 

Ftgutes  10a  and  10b  show  the  conditioned  rms  velocities,  u‘  and  v', 
plotted  against  distance  from  the  instantaneous  flame  front,  while  Table  4 
demonstrates  how  the  higher  order  moments  of  velocity  evolve  with  distance 
from  the  flame  front.  The  conditioned  velocity  statistics  were  computed  only  if 
1024  or  more  data  points  were  within  the  limia  used  to  segregate  the  data.  Note 
that,  in  general,  the  turbulent  intensities  of  both  velocity  components  show  a 
slight  decreasing  trend  as  the  flamelet  is  approached.  The  two  v-component 
intensities  at  x  •  5.8  and  6.0  cm  are  the  exception;  this  may  be  due  to  an 
interaction  with  the  wake  of  the  flame  stabilizer  rod.  This  decreasing  trend  is 
similar  to  that  reported  by  Cheng  ( 10],  though  his  measure  mens  are  conditioned 
only  on  being  within  the  reactants  or  products,  and  suggests  that  reported 
increases  in  unconditioned  turbulent  intensities  within  the  flame  brush  are  due  to 
higher  intensities  within  the  products  and  the  effects  of  measuring  intermittently 
within  the  products  or  the  reactants. 

iV  SUMMARY  AND  CONCLUSIONS 

Flamelet  position  along  a  line,  its  derivative,  and  reactant  gas  velocity 
have  been  simultaneously  measured  with  Laser  Doppler  Velocimetry  and  a 
variation  of  Laser  Tomography.  Reactant  flow  is  seeded  with  a  silicon  oil  mist 
which  provides  LDV  scattering  particles.  On  passage  through  the  flame  front 
these  particles  bum  and  evaporate:  this  property  is  exploited  to  provide  a  measure 


U-comp  turbulent  Intensity 


□  x-i.Sem 
A  x-S.Oem 
O  x-a.Sem 
x  x-y.Oem 


0-1  mm  l-tmiii  2-3mm  Wmm  4-5mm 

Distance  Irom  flamelet 

Figure  10b.  Conditioned  V-component  Turbulent  Intensities  vs. 

Distance  from  Flamelet  6  -  0.8. 

of  instantaneous  flamelet  position  within  a  measurement  volume  formed  by  the 
image  of  a  laser  beam  on  the  entrance  slit  of  a  photomultiplier  tube.  The  flamelet 
position  signal  is  differentiated,  thus  providing  information  on  the  flamelet  motion 
in  space. 

The  preliminary  results  reported  here  provide  insight  into  the  mechsnisms 
and  dynamics  of  turbulent,  V-shaped  flames.  In  particular,  it  was  found  that  the 
pdf  of  flamelet  position  is  substantially  Gaussian  and  that  the  rms  fluctuations 
increase  very  nearly  linearly  with  axial  distance.  Furthermore,  though  at  increase 
in  equivalence  ratio  increases  these  rms  fluctuations,  the  spectral  characteristics  of 
flamelet  motion  remain  unchanged.  Spectre  for  the  derivative  of  flamelet  position 
and  the  v-component  velocity  Held  are  compared,  and  demonstrate  that  between 
approximately  500  and  1500  Hz  the  flamelet  motion  is  controlled  by  the  turbulent 
velocity  field. 

Statistics  of  die  velocity  field  conditioned  on  distance  from  the  flamelet  are 
also  reported.  These  show  the  true  effect  of  the  flamelet  on  the  turbulent  velocity 
field.  It  is  found  that  fluctuating  turbulent  intensities  remain  relatively  constant  as 
the  flamelet  is  approached,  perhaps  decreeing  slightly. 

Further  application  of  this  technique  is  in  progress.  In  particular,  flamelet 
position  pdfs  at  greater  and  more  numerous  axial  locations  must  be  measured. 

Table  4a 

Moments  of  Conditional  Velocities,  +  a  0.8. 


U-taniBoncna 


Axial  Position 


U-  -,,3.8  m/scc 

Distance  from  instantaneous  flamelet  position  (mm) 


D.imm  vlmm  J.Jmn*  i  >mm  4-Smm 
Distance  from  flamelet 

Figure  10a.  Conditioned  U -component  Turbulent  Intensities  vs. 
Distance  horn  Flamelet.  ♦  -  0.8. 


5.5  cm 

0-1 

1-2 

2-3 

3-4 

4-5 

rms  (m/see) 

0.192 

0.190 

0.194 

0.198 

skew 

-0.08 

-0.01 

-0.05 

-0.04 

kurt 

3.14 

3.17 

3.14 

3.25 

6.0  cm 

rms  (m/sec) 

0.178 

0.178 

0.184 

0.186 

skew 

-0.19 

-0.08 

-0.04 

-0.07 

kurt 

3.53 

3.33 

3.21 

3.29 

6-5  cm 

rms  (m/sec) 

0.174 

0.170 

0.176 

0.181 

0.176 

skew 

-0.23 

•0.19 

•0.14 

-0.12 

-0.04 

kurt 

3.62 

3.62 

3.60 

3.36 

3.68 

7.0  cm 

rms  (m/sec) 

0.166 

0.167 

0.171 

0.174 

skew 

-0.37 

•0.33 

•0.21 

•0.16 

kurt 

3.88 

3.88 

3.45 

3-58 

*  W*\' 


Table  4b 
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V-camponmn  u_  «  nt 
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0-1 

1-2 

2-3 

3-4 

4-5 
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skew 

kurt 
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5.14 
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4.17 
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skew 
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0.16 
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kurt 
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skew 
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0.125 
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0.141 

skew 
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0.16 

0.08 

kurt 
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3.33 

aid  the  effect  of  equivalence  ratio  on  flameiet  motion  must  be  investigated  further 
Conditional  velocity  statistics  over  a  wider  range  of  axial  and  lateral  locations 
must  also  be  determined,  and  the  technique  should  be  modified  to  enable 
conditioned  velocity  measurements  within  the  products  as  welt  In  this  manner 
the  interaction  between  the  flameiet  and  the  velocity  field  can  be  fully  investigated. 
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